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TECHNICAL  SUMMARY 


In  this  report  we  describe  the  results  from  a  large  number  of 
computations  of  signal  intensity  distributions  and  channel 
symbol  bit  error  rates  for  satellite  links  under  non-Rayleigh 
fading  conditions.  The  computations  are  performed  numerically 
using  the  Fresnel-Kirchhof f  scattering  equation.  The  computa¬ 
tions  address  a  broad  range  of  parametric  representations  of 
the  ionospheric  scattering  medium;  where  applicable,  a  broad 
range  of  Fresnel  lengths  ( 4p)  is  also  considered.  The  Intent 
has  been  to  include  at  least  some  conditions  which  may  not  be 
predominant  in  the  ambient  environment  but  which  might,  per¬ 
haps,  occur  in  the  aftermath  of  one  or  more  high  altitude 
nuclear  explosions. 

The  propagation  medium  is  approximated  as  a  thin,  one-dimen¬ 
sional  (anisotropic),  spatially  structured,  phase-shifting 
screen.  Phase  screens  are  generated  as  numerical  realizations 
of  a  stationary,  spatially  correlated  Gaussian  random  process. 
The  second-order  statistics  of  the  process  are  specified  by  the 
power  spectral  density  (PSD),  which  is  the  Fourier  transform  of 
the  spatial  autocovariance  of  the  phase  screen  phase.  Thu  PSD 
is  parameterized  versus  k,  the  spatial  wavevector,  either  as  a 

I  -  n  * 

"single  power  law"  k  form  or  as  a  "double  power  law"  form 

with  an  abrupt  change  in  n  at  k  =  k_ .  An  outer  scale  size  L 

b  _  o 

is  used  to  roll  off  the  PSD  at  k  <  L  ,  and  for  n  <  3  an 

*'•  o  — 

inner  scale  size  4^  is  also  used  to  cut  off  the  spectrum  at 
large  k. 


A  k1  n  phase  screen  PSD  corresponds  to  a  k  n  three-dimensional 
refractive  index  (for  electron  density)  spatial  power  spectrum. 
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The  computed  results  are  compared  with  the  Rice  and  Nakagami-m 
models.  Under  certain  conditions  (mainly  n  <  3,5  and 
the  Rice  model  agrees  fairly  well  with  the  computed  results. 
More  generally  Rice  tends  to  be  a  worst  case.  Under  other 
specific  conditions  (mainly  n  S  4  and  with  Jl  .<<6-<<2nL  ),  the 
Nakagami-m  model  agrees  fairly  well,  but  only  at  moderate 
signal  intensity  levels;  Nakagami-m  consistently  underestimates 
the  frequency  of  deep  fades,  which  can  often  dominate  communi¬ 
cations  performance.  More  generally,  the  computed  results  show 
that  there  are  previously  unrecognized,  strong  and  systematic 
trends  in  signal  intensity  statistics  as  functions  of  the  PSD 
parameterization  and  the  Fresnel  length.  No  existing  simple 
model  will  reproduce  these  trends.  Rice  statistics  offer  an 
heuristic  "worst  case"  specification  for  generally  bounding  the 
severity  of  the  signal  intensity  effects. 

We  also  briefly  review  the  results  from  previous  comparisons  of 

ambient  environment  satellite  link  data  versus  Nakagami-m  and 

other  models.  In  these  past  studies,  initially  conflicting 

findings  have  apparently  yielded  to  a  consensus  that  Nakagami-m 

statistics  seem  to  provide  somewhat  the  better  but  imperfect 

fit  (among  those  options  considered)  to  the  ambient  environment 

data.  We  find  this  rough  consensus  to  be  consistent  with  our 

oresent  results,  since  the  conditions  (e.g.,  nS4  and  il1,<<2nL  l 

r  o 

where  Nakagami-m  best  approximates  our  computated  results  are 
also  thought  to  nominally  represent  the  most  common  features  of 
the  ambient  ionosphere.  Moreover,  the  data  used  in  these  past 
studies  did  not  accurately  sample  deep  fades,  which  are  under¬ 
estimated  by  Nakagarai-m  statistics. 
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MULTIPLY - -  BY - - -  TO  GET 


TO  GET  - - 

BY  * - 

- DIVIDE 
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**Th«  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 
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SECTION  1 
INTRODUCTION 

This  paper  addresses  signal  intensity  statistics  for  satellite 
communications  links  propagating  through  regions  of  ionospheric 
structure.  Principal  attention  is  given  to  conditions  under 
which  the  propagation  disturbances  are  not  strong  enough  to 
cause  full*  developed  Rayleigh  fading.  Also  addressed  is  the 
systematic  dependence  of  the  intensity  distribution  and  other 
signal  characteristics  upon  certain  detailed  features  of  the 
ionospheric  structure  and  tKe  Fresnel  length  ( ) . 

It  is  well-known  that  under  dynamic,  nonequilibrium  conditions 
the  ionospheric  F-region  can  become  highly  structured,  in  the 
form  of  elongated,  geomagnetic-field-aligned  "striations"  of 
excess  electron  density  (Refs.  1-4).  Transmission  of  satellite 
signals  through  such  striated  regions  can  lead  to  multiple 
scattering  and  multipath  propagation  (Refs.  5-7).  This  multi- 
path  propagation  can  cause  severe  signal  distortions,  or  "scin¬ 
tillations,"  in  the  form  of  intermittent  fading,  phase  fluctua¬ 
tions,  and  signal  angle  and  time  of  arrival  spreads  (Refs.  8-12) 
Because  of  the  massive  ionospheric  disturbances  that  would  be 
caused  by  high  altitude  nuclear  explosions,  these  effects  are 
of  particular  concern  for  military  satellite  communications 
(Refs.  13-17);  however,  they  can  also  occur  in  the  ambient 
environment  under  conditions  of  equatorial  spread-F  or  severe 
auroral  disturbances  (Refs.  18-26). 

Under  sufficiently  strong  scattering  and  multipath  conditions 
the  signal  scintillations  can  be  accurately  characterized  in 
terms  of  Rayleigh  distributed  phase  and  amplitude  (or  inten¬ 
sity)  statistics  (Refs.  27-29).  The  convenient  mathematical 
structure  of  Rayleigh  signal  statistics,  and  also  the  validity 
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of  such  statistics  over  the  broad  range  of  strong  scattering 
conditions,  have  had  important  benefits  for  the  analysis, 
design,  and  testing  of  military  satellite  links. 


Less  severe  scattering  conditions  are  also  important.  Under 
these  conditions,  the  signal  phase  fluctuations  are  typically 
rather  slow,  and  angle  or  time  of  arrival  spreads  are  generally 
negligible;  but  signal  intensity  fluctuations  are  still  signif¬ 
icant.  The  fades  are  less  deep;  but,  because  they  are  typically 
of  longer  duration,  they  can  be  very  difficult  to  mitigate.  For 
practical  applications  it  is  the  intensity  distribution  (first- 
order  statistics),  and  also  the  dynamics  (second  order  statis¬ 
tics),  of  the  intermittent  signal  fades  and  enhancements  that  is 
of  greatest  interest.  This  paper  addresses  only  first-order 
intensity  statistics  and  related  properties. 

Unfortunately,  the  proper  description  of  signal  intensity  dis¬ 
tributions  under  non-Rayleigh-fading  conditions  remains  highly 
problematical.  Several  analytical  forms  have  been  considered, 
including  the  Rice,  log-normal,  and  Nakagarai-m  distributions 
(Refs.  8,  30-33),  and  also  a  more  pragmatic  approach  to  be  dis¬ 
cussed  below.  It  seems  that  the  Rice  distribution  has  been 
favored  in  some  engineering  applications  while  the  Nakagami-m 
distribution  is  thought  to  provide  a  somewhat  better  but 
imperfect  fit  to  certain  satellite  data  taken  in  the  ambient 
environment  (Refs.  8,  31  and  Section  5). 

One  of  several  methods  which  we  will  use  for  presentation  and 
discussion  of  our  computed  results  will  be  the  channel  symbol  bit 
error  rate  characteristic  as  would  be  measured  directly  at  the 
link  demodulator  output  (see  Section  2  and  Appendix  A  for 
details).  Figure  1  depicts  predicted  channel  symbol  bit  error 
rates  (BERs)  for  differential  binary  phase  shift  key  (DBPSK) 
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versus  the  average  signal to-noise  ratio  (SNR)  for  different 
extremes  of  the  assumed  signal  intensity  distribution.  The 
region  bounded  by  the  "Slow  Rayleigh  Fading"  (SRF)  curve  por¬ 
trays  the  regime  of  Rayleigh  statistics.  The  AWGN  (additive 
white  Gaussian  noise)  curve  depicts  normal  .link  performance 
with  no  fading.  There  is  obviously  a  large  gap  in  between 
these  two  regimes.  Signal  intensity  statistics  within  this 
"gap-1  region  will  be  the  subject  of  the  present  report. 


SNR  (dB) 


Figure  1 . 


Channel  symbol  bit  error  rates  (BERs)  versus  signal- 
to-noise  ratio  (SNR  (dB))  for  DBPSK  in  Rayleigh  fading 
(SRF  curve  and  above)  and  in  non-fading  (AWGN)  regimes. 
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The  most  pragmatic  approach  for  bridging  this  gap  (Refs,  7,9,34) 

has  been  simply  to  assume  Rayleigh  statistics  whenever  the  "Rytov 

parameter"  (i.e,,  the  Rytov  approximation  to  the  signal  log- 

amplitude  variance,  herein  denoted  as  Xt  --more  fully  defined 

—y  Ry 

later)  exceeds  the  value  X‘  =0.1,  and  to  assume  negligible 

Ry_ y 

signal  distortions  whenever  X?  <  0.1.  The  rationale  for  this 

Ky 

specific  Rytov  parameter  criterion,  within  the  overall  approach, 
will  be  demonstrated  below.  An  important  practical  advantage  of 
this  approach  is  its  continuance  of  the  mathematically  convenient 
Rayleigh  specifications  into  a  broader  regime  of  application. 

One  disadvantage  under  weak  scattering  but  for  X£  >  0,1,  is  that 
it  provides  an  overiy  stressful  specification  under  ♦'hose  slow 
fading  conditions  which  are  already  quite  difficult  to  mitigate, 
and  which  also  persist  over  large  areas  and  for  long  times  in  a 
nuclear  environment.  A  further  disadvantage  is  that  it  neglects 
link  performance  degradations  altogether  when  X^  <  0.1;  these 
latter  conditions  can  occur  over  even  larger  areas,  and  for  even 
longer  times. 

Among  other  possibilities,  the  assumption  of  either  Rice  or 
Nakagami~m  signal  statistics  offers  a  means  to  more  continu¬ 
ously  bridge  the  gap  between  AWGN  (no  fading)  and  SRP  (strong 
scattering)  conditions.  Fy  definition,  both  the  Rice  and 
Nakagami-m  models  are  parameterized  by  the  scintillation  index 
S^,  which  is  the  signal  intensity  variance.  As  ■*  1.0,  both 
Rice  and  Nakagami-m  approach  the  SRP  limit.  As  s|  0,  both 
approach  the  AWGN,  or  no-fading  limit.  In  between,  however, 
the  two  models  may  differ  appreciably. 

Por  weak  scattering,  3  4X*  .  Included  in  Figure  2  are  the 

BER  vs.  SNR  characteristics  (again  DBPSK )  predicted  from  both 

Rice  and  Nakagami-m  statistics,  and  specifically  at  S*  =  0.378, 

which  also  corresponds  to  xf  =0.1  for  a  particular  scattering 

Ky 
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medium  parameterization  of  common  interest  (see  Table  1  of  Sec¬ 
tion  3).  It  can  be  seen  that  either  curve  lies  approximately 
midway  between  the  SRF  and  AWGN  extremes.  Observations  of  this 

nature  provide  the  rationale  for  selecting  X*  *  0.1  as  a  conven- 

xy 

ient  dividing  line  between  Rayleigh  fading  and  no  fading  in  the 
simpler  pragmatic  approach  described  above.  However,  it  can  also 
be  seen  from  Figure  2  that  even  the  BER  vs.  SNR  characteristics 
of  Rice  versus  Nakagami-m  statistics  may  differ  appreciably  under 
conditions  of  practical  interest.  As  described  below,  comparable 


SNR  (dB) 

Figure  2.  Slow  Rayleigh  fading  (SRF)  and  non-fading  { AWGN ) 
BER  vs.  SNR  characteristics  (from  Figure  1),  as 
compared  with  Rice  (R)  and  Nakaganii-m  (N)  statis¬ 
tics  at  s|  =  0.378  --  all  for  DBPSK . 
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or  even  larger  differences  can  also  occur  when  the  models  are 
compared  with  more  detailed  theoretical  calculations.  These 
differences  further  help  to  motivate  the  present  investigation. 

In  this  report,  we  will  describe  the  results  of  a  theoretical 
and  computational  investigation  of  signal  .intensity  statistics 
under  non-Rayleigh  fading  conditions.  A  large  number  of  compu¬ 
tations  of  received  signal  charactertics  have  been  performed, 
using  the  Fresnel-Kirchhof f  scattering  equation,  and  encompass¬ 
ing  a  fairly  broad  range  of  parametric  descriptions  of  the 
striated  ionospheric  propagation  medium.  Strong  and  systematic 
trends  can  be  identified  in  the  results,  as  summarized  below. 

We  find  chat  neither  Nakagami-m  nor  Rice  statistics  are  reliable 
in  general.  Within  the  parametric  regimes  of  greatest  present 
interest,  the  Rice  distribution  tends  to  provide  an  approximate 
upper  bound  on  the  frequency  of  deep  fades,  and  on  predicted 
channel  symbol  BERs ,  while  the  detailed  BER  vs.  SNR  curves  may  be 
either  higher  or  lower  than  predicted  by  Nakagami-m  statistics, 
depending  both  on  the  parametric  model  of  the  scattering  medium 
and  on  the  problem  SNR  and  Fresnel  length. 

We  consider  situations  in  which  the  ionospheric  electron  density 
or  refractive  index  structure  can  be  parameterized  in  terms  of  a 
three-dimensional  spatial  power  spectrum  approximately  of  the  form 
k~n,  where  k  -  2n/length  is  the  spatial  wavevector.  By  defini¬ 
tion,  s  =  n- 2  is  then  the  "spectral  index,"  and  in  the  thin  phase 
screen  approximation  (Section  2)  the  phase  screen  PSD  varies  as 
k1_n.  In  the  chosen  parameter izat ions ,  the  exponent  n  may  be 
constant  over  a  broad  range  of  wave  vectors,  or  it  may  be  allowed 
to  change  discontinuous!'/  at  k  =  kg  (i.e.,  a  k  form  at  k  <  kg, 
and  a  k  n2  form  at  k  >  kg) .  in  either  case,  we  find  that  the 
signal  intensity  statistics,  BER  vs.  SNR  characteristics,  and 
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o^hei  results  as  well  can  be  at  least  approximately  characterized 
.n  terms  of  an  "effective"  value  of  n  which  provides  a  reasonable 
k-‘1  fit  to  the  spectrum  for  wave  vectors  generally  somewhat  less 
than,  and  in  the  vicinity  of,  k  £  2rr/flp,  a  more  precise  charac¬ 
terization  is  provided  by  the  specific  examples  given  later. 

The  findings  are  partially  illustrated  in  Figure  3,  which  pre¬ 
sents  calculated  BER  vs.  SNR  characteristics  (again  for  DBPSK)  as 
a  function  of  the  exponent  (n)  of  an  approximately  k1"^  phase 


SNR  (dB) 


Figure  3.  BER  vs.  SNR  characteristics  for  DBPSK,  as  calcu¬ 
lated  from  a  k_n  refractive  index  PSD  with  n  =  3.0, 
3.5,  4.0,  4.5  and  5.0 
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screen  PSD  all  still  at  =  0.378.  Also  shown  are  the  corres¬ 
ponding  SRF  and  AWGN  limits.  For  n  £  3.5,  the  resulting  BER  vs. 
SNR  characteristics  fairly  clcseiy  approximate  tho^e  nredic  ed  by 
the  Rice  intensity  distribution  (compare  Fig.  2).  For  n  3  4.0 
the  BER  versus  SNR  characteristics  lie  intermediate  between  those 
from  Rice  and  Nakagami-m  statistics,  favoring  Nakagami-m  only 
at  the  lower  SNR  levels.  For  n  ;>  4.5,  the  BER  characteristics 
lie  below  those  from  Nakagami-m  statistics  at  lower  SNR  levels, 
but  above  Nakagami-m  at  higher  SNR  levels. 

Related  strong  and  systematic  trends  are  found  in  the  signal 
intensity  distributions  and  in  the  versus  S£  relationships 

as  functions  of  the  parameterization  of  the  scattering  medium  and 
the  Fresnel  length,  both  for  single  and  double  power  law  phase 
screen  PSDs .  In  addition  to  their  relevance  for  the  design  and 
analysis  of  military  satellite  communications  links,  some  of 
these  results  may  also  be  applicable  to  the  interpretation  of 
satellite  data  taken  in  the  ambient  environment. 
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SECTION  2 

MATHEMATICAL  BACKGROUND  AND  APPROACH 


In  this  section  we  describe  the  mathematical  background  and 
general  parameterization  schemes  of  interest  for  the  present 
investigation.  This  will  include:  the  mathematical  descrip¬ 
tion  of  the  ionospheric  scattering  medium  (subsection  2.1);  the 
numerical  generation  of  specific  realizations  of  the  medium 
(subsection  2.2);  the  Fresnel -Ki rchhof f  formalism  for  calcula¬ 
ting  the  properties  of  received  signals  (subsection  2.3);  and 
the  definition  of  various  measures  of  signal  scattering  intens¬ 
ity  and  other  key  parameters  such  as  signal  line-of-sight.  (LOS) 
phase  variance  (or|),  the  scintillation  index  (S^),  the  Rytov 
parameter  approximation  to  the  log-amplitude  variance  (X^  ) , 
and  the  Fresnel  length  -~a.il  in  subsection  2.4.  In  subsec¬ 

tion  2.5  we  describe  the  general  properties  of  Rice,  Nakagami- 
m,  and  Rayleigh  intensity  statistics.  In  subsection  2.6  we 
describe  the  general  approach  for  representing  slow  fading 
effects  in  terms  of  BER  vs.  SNR  characteristics  for  typical 
digital  modems,  with  SERs  as  measured  at  the  basic  channel 
symbol  level,  without  the  benefits  of  error  detection-cor¬ 
rection  coding  and  interleaving.  Additional  details  are  also 
provided  in  Appendixes  A  and  B. 

2.1  THE  SCATTERING  MEDIUM. 

In  the  present  investigation,  the  striated  ionospheric  scat¬ 
tering  medium  will  be  represented  as  a  spatially  thin,  phase- 
shifting  screen.  The  partial  neglect  of  certain  finite-medium- 
depth  effects  on  signal  propagation  is  an  approximation,  of 
course,  but  has  been  assessed  for  its  reliability  in  Reference 
34  Moreover,  the  overall  findings  from  the  subsequent  analy¬ 
sis  are  clearly  sufficiently  robust  as  to  be  qualitatively 
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unaffected  by  this  approximation.  The  signal  phase  shift 
(more  precisely,  its  deviation  relative  to  the  mean)  is  regar¬ 
ded  as  a  Gaussian  random  variable  in  the  (x,  y)  coordinates  of 
the  plane  perpendicular  to  the  signal  1 ine-of-sight  (LOS).  The 
phase  statistics  are  therefore  fully  defined  by  the  autocovariance, 

<4>(x,  y)<t>(0,0)>  =  B  ( x  ,  y )  .  (1) 

We  effectively  neglect  the  mean  phase  shift,  phase  advance,  and 
group  delay  by  setting  <^>  »  0.  In  practice,  we  will  also 
neglect  the  y-dependence  of  £(x,y)  and  B(x,y).  This  last 
assumption  is  expected  to  have  very  little  qualitative  effect 
on  the  signal  scattering  statistics;  for  a  related  discussion, 
see  Reference  34. 

The  medium  is  specifically  parameterized  in  terms  of  the  phase 
screen  power  spectral  density  (PSD),  defined  as 

4>(k)  =  f  B(x)  exp(-ikx)  dx.  (2) 


Several  different  PSD  parameterizations  are  addressed.  The 
simplest  is  a  "single  power  law"  of  the  form  (for  v>l) 


♦  (k) 


2  /n  » 


H  v) 


*  r(v-*)  [1+Llk2)v 


(3) 


2 

Here,  er£  is  the  phase  screen  phase  variance  (as  sampled  over  an 
ensemble  of  LOS  ray  paths),  and  LQ  is  an  "outer  scale  length" 
effectively  defining  the  maximum  striation  scale  size  as 
measured  perpendicular  to  the  LOS.  The  exponent  v  is  related 
to  the  more  conventional  spectral  index  as  v  =  (s+l)/2.  This 
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parameterization  of  the  phase  screen  PSD,  essentially  as  a  k~2l/ 
power  law  (i.e,,  for  k>>L0-1)’  ^lso  implies  a  parameterization 
of  the  underlying  striation  electron,  density  three-dimensional 
PSD  as  a  k'n  form,  with 

n=s+2=2v+l.  (4) 

The  special  case  v  =  1  (also  s  =  1,  n  =  3)  requires  introduction 
of  an  additional  "inner  scale  length"  cutoff  (.0^)  so  that  the 
underlying  electron  density  variance  will  be  finite.  Therefore, 
for  v  -  1  we  use 


+  (k) 


2  ^  2  L 

- Z  KX(Z)  exp^i/L,.,) 

1  +•  L^k^ 


(5) 


where  Z2L2  =  02(1+L2k2),  and  K1(Z)  is  a  Bessel  function. 


We  will 
tion  of 


also  address  a  more  general  "two  power  law" 
the  form 

/ 

.2 

,  k  <  kB, 


Lo** 


4>(k)  = 


(1  +  hlv.2)  1 


( L2k2 ) V  ^ 


,  k  >  kB. 


parameteriza- 


(6) 


Conditions  of  continuity  and  normalization  determine  Cj  and  C2 
as  functions  of  v^,  v2,  and  LQkB .  These  relations  are  detailed 
in  Appendix  A.  Evidently,  k0  defines  a  length  scale  (2n/k0)  at 
which  there  is  a  "break"  in  the  spectral  index.  Under  condi- 


tions  cf  .Interest  LQ  ~  10  km  cr  so,  while  2n/kg  is  typically  on 
the  order  of  several  hundreds  of  meters  (Refs.  35-36).  Thus, 
when  =  v2  the  two-power  law  PSD  very  closely  approximates 
the  one-power  law  form.  As  for  the  spectral  index,  many  condi¬ 
tions  of  practical  interest  are  presently  thought  to  be  encom¬ 
passed  by  (Ref.  37) 

1  *  i  3/2  *  *2  *  $ / 2 .  (7) 

Even  more  complex  "three  power  law"  models  are  also  under 
present  consideration,  based  on  recently  acquired  ambient  iono¬ 
sphere  satellite  data  (Ref.  37);  this  report  will  emphasize 
effects  which  are  largely  invariant  to  such  further  details. 

2.2  NUMERICAL  REALIZATIONS. 

Given  these  descriptions  of  the  scattering  medium  in  terms  of  a 
thin  phase  screen,  with  the  LOS  phase  treated  as  a  spatially 
stationary  Gaussian  random  process  having  a  specified  PSD, 
representative  phase  screen  realizations  can  then  be  numeric¬ 
ally  generated  by  standard  statistical  sampling  methods.  A 
representative  sampling  of  such  realizations  can  then  be  used 
to  numerically  calculate  the  statistical  properties  of  the 
received  signal . 

For  each  defined  phase  screen  PSD  parameterization,  a  set  of 
specific  phase  screen  realizatons  can  be  generated  as 


N-l 

$(xm)  =  Re  ^  Dnrn  exp {  2ninm/N )  j  (8a) 

n=0 
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where 


2n/NAx 


2n/xN, 


(8b) 


Ak 

Dq  =  +  ( k=0 )/2n, 
D^o  =  +(nAk)/n. 


(8c) 

(8d) 


The  quantity  rn  =  rin  +  ir2n  is  u  c<r>mPiex  Gaussian  random  vari¬ 
able  with  <rx>  =  <r2>  -•  0,  <r\>  =  <r|>  =  1,  and  <r1r2>  =  0. 

The  numerical  and  implementat ional  details  are  further  defined 
in  Appendix  B . 


2.3  FRESNEL-KIRCHHOFF  SCATTERING  FORMALISM. 


With  neglect  of  large  angle  scattering,  geometric  divergence, 
absorption,  and  antenna  gain  effects,  the  received  signal  can 
be  determined  from  the  Fresriel-Kirchhof  f  scattering  equation 
as  (Refs.  34,38) 


h  ( x ) 


(\Z 


-1/2 


-  i  n  /  4  r°° 

e  j  dx ' 


—  CO 


fci*(x' ) exp( 2ni ( x-x ' )2/2XZ*]  . 

(9) 


Here,  h(x)  is  a  complex  received  signal  amplitude  modulation  a 
position  x  in  the  receiver  plane  ( y-dependent  variations 
neglected).  As  usual,  X  is  the  signal  wavelength.  The  func¬ 
tion  <f>(x' )  is  a  specific  realization  of  the  phase-shifting 
screen,  determined  as  defined  above.  The  parameter  Z*  is  an 
equivalent  signal  path  length  defined  as 

ZRZT 


Z* 


+  Z  * 


(  10) 


where  ZR  is  the  distance  to  the  receiver  from  the  center  of  the 
phase-shifting  medium,  and  ZT  is  the  corresponding  transmitter- 
to-medium  distance.  Notice  the  importance  of  the  parameter 
combination  \Z* ,  which  is  one-half  the  square  of  the  Fresnel 
length  (J}R).  The  role  of  the  Fresnel  length  will  further 
discussed  in  subsection  2.4. 

Given  a  specific  phase  screen  realization,  as  <P(Xj)  =  Ax), 
we  can  numerically  evaluate  Equation  9  to  determine  the  x-depen- 
dence  of  the  received  complex  signal  amplitude  h(x).  From  suffi¬ 
ciently  numerous  calculations  of  this  sort,  and  including  a 
sufficiently  numerous  set  of  phase  screen  realizations,  the  sta¬ 
tistics  of  the  received  signal  intensity  I  =  |h|2  are  ultimately 
determined.  The  numerical  version  of  Equation  9  Is  written  as 


h(  x, 


AX 

VTT*- 


-  in/4 


J+N/2-1 

5  • 

m=J-N/2 


14>(x 


m' 


exp  ^2ni  (  Ax)  2  (  J>-m  ) 2/ 2XZ  * j 

(11) 


Additional  details  are  given  in  Appendixes  A  and  B. 

2.4  OTHER  KEY  PARAMETERS. 


The  total  strength  of  the  scattering  medium  can  be  charac- 
terized  in  terms  of  «r£,  which  is  the  signal  phase  variance  as 
would  be  measured  over  an  ensemble  of  LOS  ray  paths  fully 
sampling  the  scattering  medium.  Under  certain  circumstances, 
the  LOS  signal  phase  variance,  which  normalizes  the  thin  phase 
screen  PSD  (i.e.,  at  k  -  0),  can  be  related  to  the  non-thin 
striated  medium's  electron  density  variance  as 


cr 


2  X. 2 


r(v~i/2) 
n  v-i ) 


(  12  ) 
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where  or2  is  the  electron  density  variance,  and  L  is  the 
effective  thickness  of  the  scattering  medium.  LQ  is  now  the 
effective  outer  scale  length  specifically  as  measured  along 
the  signal  LOS;  re  =  2.82  x  10_1®  m  is  the  classical  electron 
radius.  This  simple  relationship,  which  pertains  only  to  the 
specia.  case  of  a  single-power  law  k'2l/  PSD  (with  v  >  1), 
nonetheless  helps  to  define  the  general  nature  of  the  relation¬ 
ship  between  the  LOS  signal  variance  of  a  thin  phase  screen 
model  and  the  electron  aensity  variance  of  a  f ini te- thickness 
striated  ionospheric  medium.  Additional  details  are  given  in 
Appendix  A. 

The  intensity  of  signal  scattering  effects  can  already  be  seen 
from  Equation  9  to  depend  not  only  upon  cr|  but  also  upon  the 
correlation  of  these  phase  fluctuations  over  length  scales 
measured  relative  to  the  wavelength-geometry-dependent  Fresnel 
length.  Because  of  this  complication,  other  single-parameter 
measures  of  the  strength  of  the  signal  perturbation  effects 
have  traditionally  been  introduced,  as  defined  below. 


The  signal  scintillation  index,  s| ,  is  defined  as  the  variance 
of  the  signal  intensity 


S 


2 

4 


<I2> 


-  <I> 


<I> 


(  13  ) 


Recall  that  the  intensity  is  determined  as  I  =  |h|2,  where  h(x) 
is  the  complex  signal  amplitude,  as  given  in  the  Fresnel- 
Kirchhoff  equation.  Obviously,  for  no  scattering  <I‘>  =  <I>  , 
and  =  0.  Under  fully  developed  strong  scattering  condi- 
tions,  Rayleigh  fading  is  rapidly  approached,  with  £  1.0 
(Refs.  27-29).  Under  less  than  fully  saturated  Rayleigh  fading 
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conditions,  the  scintillation  index  can  also  temporarily  "over¬ 
shoot1'  (i.e,,  s|  >  1-0),  to  a  degree  which  depends  upon  the 
specific  parameterization  of  the  scattering  medium,  and  which 
also  will  be  detailed  by  the  present  investigation. 

A  second  commonly  used  measure  of  the  strength  of  signal  scat¬ 
tering  is  the  Rytov  approximation  to  the  log-amplitude  vari¬ 
ance,  herein  defined  as 

%  “  ♦OOsanZlZ.^M,,).  ,14) 

For  weak  scattering  conditions  it  can  be  shown  that  S  4 

(see  Appendix  A).  However,  for  increasingly  strong  scattering 

conditions  X„  increases  without  limit  and  ceases  to  be  a 
Ky 

physically  meaningfully  approximation  to  the  true  log-amplitude 
variance,  while  S2  ultimately  equilibrates  at  =  l.o.  The 
detailed  relationships  between  and  will  be  further  dis¬ 

cussed  in  the  following  sections. 

Both  Equations  9  and  14  demonstrate  the  critical  role  of  the 
signal-geometry-dependent  Fresnel  length,  defined  as 

*  \  /  2 

iF  -  [2VZ  r  .  (15) 

The  importance  of  the  Fresnel  length  in  these  problems  can  be 
more  fully  appreciated  by  rewriting  Equation  9  as 


h(y*F) 


-  in/4 
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where  y  =  x/JJp  and  y'  =  x'/JJp.  Thus,  the  first  order  received 
signal  statistics  depend  only  upon  the  total  strength  of  the 
phase  screen  fluctuations,  as  measured  by  cr^,  and  upon  their 
variations  upon  length  scales  measured  in  terms  of  ,#F. 

2.5  SIGNAL  INTENSITY  STATISTICS  MODELS. 

In  lieu  of  detailed  calculations  of  the  types  to  be  described 
in  the  present  report,  several  different  models  have  been  pro¬ 
posed  for  signal  intensity  statistics  under  conditions  of  iono¬ 
spheric  scattering  and  multipath.  Most  extensively  studied,  by 
far,  has  been  the  Rayleigh  distribution.  It  has  been  estab¬ 
lished  that  under  sufficiently  strong  scattering  conditions  the 
signal  intensity  probability  distribution  will  approach  the 
Rayleigh  limit,  defined  as  (Ref.  39) 

P(I)  =  e-1  (17) 

where  I  is  the  (normalized)  signal  intensify,  and  P(I)  its  prob¬ 
ability  distribution. 

The  Rice  distribution  assumes  implicitly  that  the  signal  can  be 
decomposed  into  the  sum  of  a  Rayleigh  component  and  an  unscat¬ 
tered  component.  Generalization  to  Include  superimposed  uncor¬ 
related  phase  effects  is  also  fairly  straightforward  for  both 
Rayleigh  and  Rice  statistics.  The  Rice  distribution  is  defined 
(Ref.  39)  by 


with 


Pd)  =  «e~*(I!'1)+1  i0(2v«(«-rrr  ) 
«=(!-(!-  s| ) 1 /2  ]  " 1 . 


(18) 
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When  3 1  -*  o,  oc  -*  eo,  and  L5 ( I )  -*  S(I  --  1);  with  $(x)  being  the 
Dirac  delta  function.  When  s|  *♦  1 ,  *  -*  1  ,  and  P(I)  approaches 
the  Rayleigh  limit.  Notice  that  for  s|  >  1.0  the  entire  pre¬ 
scription  of  Rice  (and  also  Rayleigh)  signal  statistics  is 
meaningless.  Plausible  conditions  leading  to  s|  >  1.0  will  be 
described  in  this  report;  in  practice  however,  these  conditions 
do  not  appear  to  be  appreciably  more  stressing  for  communica¬ 
tions  links  than  either  Rice  or  Rayleigh  statistics. 


The  Nakagami-m  distribution  is  defined  in  terms  of  m 


as 


P(  I) 


r(m) 


(  IS) 


Again,  when  s|  -*  o,  r<»  -»  <*>;  and  the  distribution  converges  to 
£(I-1).  When  =  1.0,  m  =  1.0;  and  Rayleigh  statistics  are 
regained.  Unlike  the  Rice  distribution,  there  seems  to  be  no 
underlying  physical  model,  valid  or  otherwise,  for  Nakagami-m 
statistics.  Although  approximate  plausibility  arguments  have 
been  presented  (Ref.  40),  the  model  is  basically  empirical;  and, 
as  we  shall  see,  its  erstwhile  "success"  in  approximately  fitting 
ambient  environment  satellite  data  seems  to  have  been,  at  least 
in  part,  a  coincidence  of  nature  and  of  the  limitations  on  the 
candidate  models  and  data  sets  employed. 

A  further  principal  difference  between  Rice  and  Nakagami-m  is 
that  the  Rice  P(T)  versus  1  distribution  is  always  finite  at  I 
=  C  (except,  of  course,  in  the  -»  o  limit);  whereas,  the 
Nakagami-m  distribution  always  vanishes  at  I  =  0  (except  in  the 
limit  when  s|  *♦  1.0).  Thus,  Rice  generally  gives  a  higher 
probability  of  very  deep  fades  than  does  Nakagami-m.  In  the 
opposite  extreme,  Rice  also  generally  gives  a  higher  probabil¬ 
ity  of  strong  signal  enhancements  (I  >>  1.0)  than  Nakagatni-ra. 
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Because  of  the  high  instantaneous  BERs  encountered  during  deep 
fades,  the  Rice  intensity  distribution  provides  the  more  stress¬ 
ing  specification  of  the  two. 

2.6  ERROR  RATE  CHARACTERISTICS. 

Given  a  determination  of  the  signal  intensity  distribution 
P(I),  whether  based  on  one  of  the  above  models  or  on  the 
detailed  Presnel-Kirchhof f  solutions,  it  is  then  straight¬ 
forward  to  determine  the  channel-symbol-level  BER  versus  SNR 
characteristics  under  the  assumption  that  the  fades  (and  the 
related  signal  phase  distortions)  are  negligibly  slow  compared 
to  the  link  modulation,  and  also  as  compared  to  the  response 
rates  of  any  receiver  tracking  loops. 

This  slow-fading  channel  symbol  BER  is  simply  the  long-term 
average  of  the  instantaneous  demodulator  output  BER  (i.e,, 
witnout  error  correcting  coding)  versus  the  link's  instantan¬ 
eous,  fluctuating  SNR  characteristics  as  averaged  over  the 
actual  distribution  of  instantaneous  signal  intensities,  and 
including  both  fades  and  also  intermittent  signal  enhancements, 

<Pe(Y)>  *  f  d I  P ( I )  Pe ( Y I ) •  (20) 

*  o 

Here  <Pe(Y)>  is  the  average  BER,  and  y  is  the  average  SNR. 

P(I)  is  the  signal  intensity  distribution  for  the  fading  chan¬ 
nel,  and  Pe ( yl )  is  the  unperturbed  BER  versus  SNR  characteris¬ 
tic  (i..e.,  for  AWGN)  at  an  instantaneous  SNR  -  yl. 

In.  this  report,  BER  versus  SNR  characteristics  are  computed  as 
a  function  of  the  specific  scatt*  ring  medium  and  Fresnel  length 
parameterizations  and  for  a  variety  of  common  binary  digital 
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modems.  These  include:  coherent  phase  shift  key  (CPSK),  dif¬ 
ferentially  encoded  but  coherently  demodulated  phase  shift  key 
(APSK),  differentially  encoded  and  demodulated  phase  shift  key 
(DPSK),  and  frequency  shift  key  (FSK).  Results  are  also  calcu¬ 
lated  for  quaternary  and  8-ary  FSK.  The  qualitative  results 
are  understandably  somewhat  similar  in  all  cases;  therefore,  in 
the  main  text,  the  special  case  of  DBPSK  will  be  used  to  illus¬ 
trate  the  trends.  A  further  discussion  will  be  found  in  Appen¬ 
dix  A,  and  a  compilation  concerning  other  modems  will  also  be 
provided  in  Appendixes  E  and  I.  For  now,  we  note  that,  for 
DBPSK, 

Pe(rD  “  1/2  exp(-Tl) •  (21) 
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SECTION  3 

SINGLE  POWER  LAW  RESULTS 


The  forms  for  the  single  power  law  PSD  that  we  have  employed  in 
our  calculations  are  given  in  Section  2  by  Equation  3  for  v  >  1 
and  Equation  5  for  v  =  1.  As  already  suggested,  for  2n  LQ  >> 

>>  (where  A ^  2  0  unless  u  =  1)  the  results  for  single  power 
laws  are  quite  insensitive  to  the  value  of  the  Fresnel  length  and 
essentially  depend  only  on  the  PSD  exponent  and  the  strength  of 
the  scattering.  The  bulk  of  calculations  were  performed  for  a 
Fresnel  length  of  Jlp  =  693  m,  and  LQ  =  10  km.  Twelve  values 
for  were  selected  in  the  range 

0.01  $  X|  $  20.0. 

These  selected  xf  values  essentially  span  the  range  from  negli- 

Ky 

gible  fading  up  through  saturated  Rayleigh  fading. 

Here  we  will  discuss  the  general  behavior  of  the  S4  -versus- 
Xpy  relationship,  the  intensity  probability  distributions,  and 
the  channel  symbol  bit  error  rates  that  result  from  single 
power  law  striation  PSDs .  The  behavior  will  be  illustrated 
with  selected  examples;  numerous  additional  calculated  results 
can  be  found  in  the  Appendixes.  Appendix  C  contains  the  inten¬ 
sity  probability  distributions  compared  either  to  Rice  and 
Nakagami-m  (for  s|  <  1 )  or  to  Rayleigh  (for  S4  ?  1).  Appendix  D 
contains  deep  fade  behavior  for  selected  distributions  compared 
to  Rice  and  Nakagami-m.  Appendix  E  presents  BER  curves  for  the 
selected  modems  for  0.01  $  Xpy  <  0.4.  Finally,  the  comparison 
of  the  calculated  BFRs  for  DBPSK  to  the  corresponding  Rice  and 
Nakagami-ra  behavior  is  given  in  Appendix  F. 
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3.1  SIGNAL  CHARACTERISTICS. 


3.1.1  Scintillation  Index. 

The  scintillation  index  ,  Equation  13,  has  been  calculated  for 
all  the  selected  values  of  5wj|  and  s  *  n--2  =  2v-l .  These  results 
are  giveu  in  Table  1  and  presented  graphically  in  Figure  4.  Some 
well-known  behavior  is  evident.  For  small  values  of  the 

results  '-'or  all  PSDc  satisfy  the  approximate  linear  relationship 
Gt  s  4-X^  .  For  s  =  1,  S4  mono  topically  increases  and  saturates 
at  unity.  For  somewhat  larger  values  of  s,  there  is  a  modest 
overshoot  e.nd  then  a  relaxation  down  to  unity.  For  still  larger 
values  of  s  (s  >  2),  there  can  be  a  significant  overshoot  of 
unity,  with  the  relaxation  back  to  unity  occurring  only  at  quite 
large  valu.es  of  xjj-y .  This  behavior  is  due  to  the  relatively 
coherent  signal  focusing  that  can  be  caused  by  long  wavelength 
ionospheric  disturbances,  which  are  statistically  more  prominent 
for  steeper  PSD  distributions. 


Table  i.  s|  vs  and  s. 


CASE 

*Ry 

s  =  1 

s  -  1.5 

s  =  2 

S  »  2 . 5 

s  =  3 

a 

0.01 

0 . 039 

0 . 039 

0.040 

0.040 

0.041 

b 

0.025 

0.095 

0 . 037 

0.099 

0.103 

0. 109 

c‘ 

0.05 

0 . 1  ol 

0  186 

0.197 

0.213 

0 . 24, 

d 

0. 1 

0 . 327 

0 . 343 

0.378 

0.438 

0 . 548 

e 

0.25 

0 .616 

0 . 664 

0 . 764 

0.959 

1  .  28 

f 

0 . 4 

0.763 

0.825 

0.962 

1 . 21 

1 . 63 

9 

0.7 

0 . 899 

0 . 958 

1 .12 

1.42 

1 . 96 

h 

1 . 5 

0.975 

1 . 03 

1 . 20 

!  .  55 

2  .  15 

j. 

3 . 0 

0 . 992 

1  .  03 

1  .  20 

1 . 57 

2  13 

j 

7 . 0 

0 . 997 

1 . 02 

1 .17 

1 . 53 

1 . 99 

k 

10 . 0 

0 . 999 

1 . 02 

1  .  16 

1 . 49 

1 .90 

1 

20.0 

0 . 998 

1 . 01 

1 .13 

1.41 

1.71 
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„<ytov  Parameter  (  XRy> 


Figure  4.  Scintillation  Index  (s|)  versus  the  Rytov  parameter 
— ^ 

( XRy )  for  various  single  power  law  PSDs. 

There  has  been  some  confusion  in  the  research  literature  as  to 
the  behavior  of  at  values  greater  than  unity.  Thus,  for 
example,  Rlno  and  Owen  (Ref.  10)  state  that  for  s  >  2  and  in 
the  strong  scatering  limit  will  converge  to  a  greater  than 
unity  value  of  (6-s)/(4-s).  However,  the  evidence  to  this 
effect  presented  by  Reference  10  and  others  cited  therein  has 
not  seemed  persuasive;  and  our  present  results  clearly  show 


that  the  purported  behavior  does  not  occur.  We  have  not  car¬ 
ried  the  computations  to  sufficiently  strong  scattering  to 

O 

verify  that  Sj  always  converges  exactly  to  one,  but  we  also  can 
offer  no  reason  to  suppose  otherwise. 

The  results  in  Table  1  and  Figure  4  reveal  a  strong  and  system- 
matic  trend  in  the  s|  -versus-  relationship  as  a  function 

of  the  single  power  law  PSD  spectral  index.  There  is  an  even 
stronger  and  more  basic  trend  which  should  also  be  noted, 
having  to  do  with  the  effect  of  the  spectral  index  on  the 
relationship  of  s|  to  the  standard  deviation  of  the  scattering 
medium's  in-situ  electron  density  or  (equivalently)  refractive 
index  fluctuations.  Unlike  Xpy,  which  is  a  theoretically 
derived  quantity,  studies  of  the  performance  of  military  satel¬ 
lite  links  in  nuclear-perturbed  environments  typically  use  the 
electron  density  standard  deviation  aNe  to  characterize  the 
basic  environment.  The  electron  density  spatial  power  spec¬ 
trum,  which  then  essentially  determines  the  PSD  in  the  thin 
phase  screen  approximation,  is  separately  specified  as  an  over¬ 
lay  to  the  basic  environment,  with  the  assumed  PSD  details 
being  subject  to  a  range  of  parameterization  uncertainties  such 
as  those  addressed  in  this  report. 

For  s  >  1  the  relationship  between  X^y  and  is  formally 

defined  by  Equations  12  and  14  of  Section  2.  For  the  LQ  and  ip 
values  used  here,  the  X^y  -  to  -  proportionality  has  been 
numerically  evaluated  and  is  given  in  Table  B--1  of  Appendix  B. 
Using  these  results  together  with  those  of  Table  1  or  Figure  4 
the  corresponding  s|  -versus--  erNe  relationship  can  be  found. 

The  result  is  depicted  in  Figure  5.  There,  «rNe  is  given  in 
relative  units,  parametric  in  L0 ,  LQ '  ,  ip,  I,,  and  X.  . 
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It  can  be  seen  from  these  calculated  results  that  the  range  of 
»Ne  values  needed  to  drive  the  weak- to-strong  scattering  tran¬ 
sition  can  be  a  very  sensitive  function  of  the  PSD  spectral 
index.  The  same,  of  course,  is  equally  true  of  the  regime  of 
validity  of  strong  scattering  and  Rayleigh  signal  statistics. 
This  sensitivity  to  the  assumed  PSD  spectral  index  can  be  a 
very  important  consideration  of  assessments  of  military  satel¬ 
lite  link  performance  in  nuclear  -  perturbed  environments. 


Figure  5.  -versus-  *Ne  relationship  ( orNe  relative)  for 

single  power  law  PSDs  with  different  spectral  indexes. 


With  double  power  law  PSD  parameterlzatlons  (as  discussed  in 
Section  4)  these  dependences  will  be  more  complex,  and  the 
Fresnel  length  will  also  frequently  be  an  important  parameter. 
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Moreover,  the  sensitivity  to  either  one  of  the  two  spectral 
indexes  of  a  two  power  law  PSD  (i.e.,  with  the  other  fixed) 
will  typically  be  somewhat  less  than  portrayed  by  Figure  5. 

3.1.2  S4  <  1  Intensity  Distributions. 

We  will  now  present  a  representative  set  of  calculated  results 
for  signal  intensity  distributions  (i.e.,  P(I)  versus  I)  and 
compare  these  results  with  the  Rice  and  Nakagami-m  models. 
Emphasis  here  will  be  on  behavior  for  S4  <  1.  Specifically, 
we  will  look  at  X^y  =  0.1,  which  corresponds  to  s|  in  the 
range  from  0.327  to  0.548  (see  Table  1),  depending  on  the 
value  of  the  PSD  exponent.  This  choice  of  lies  inter¬ 

mediate  between  the  AWGN  and  SRF  conditions,  and  therefore 
helps  to  bring  out  more  clearly  the  various  detailed  differ¬ 
ences  between  the  models  and  computations.  Of  course,  the 
model  predictions  change  somewhat  as  s|  changes,  even  though 
X^y  is  fixed.  Numerous  additional  results  at  other  Xj|y  values 
(but  still  for  S4  <  1 )  can  be  found  in  Appendixes  C  and  D,  with 
Appendix  D  emphasizing  the  behavior  at  deep  fades. 

Figure  6  shows  the  computed  and  model  predictions  at  X^y  =  o.l 
for  an  s  =  1  (also  n  =  3,  v  =  1 )  PSD.  Figure  6a  gives  the 
results  at  moderate  intensity  values,  and  Figure  6b  gives  the 
results  for  deep  fades.  Throughout  this  report  the  average  (or 
unperturbed)  signal  intensity  is  <I>  =1.  It  is  immediately 
apparent  that  Rice  statistics  provide  an  excellent  fit  to  the 
computed  results  in  this  case  (s  =  1).  This  is  rather  as 
expected,  since  an  s  =  1  PSD  gives  greater  weight  to  shorter 
spatial  wavelength  ionospheric  structures,  which  thus  tends  to 
put  the  receiver  in  the  far  zone  from  the  scatterers.  Under 
far  zone  conditions.  Rice  statistics  «ire  generally  expected  to 
be  mcst  accurate  (Ref.  41). 
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Figure  6.  Comparison  of  computed  (C)  results  with  Rice  (R)  and 
Nalcagaml-m  (N)  models  at  »  -  i,  S*y  .  0>1<  s2  .  0_3„ 

By  comparison,  we  see  that  Nakagaml-m  provides  a  poor  fit  to 
the  computed  results  under  these  conditions.  It  gives  tcc  high 
a  probability  of  near-nominal  (l.e.,  I  a  1.0)  signal  Intensi- 
ties,  and  too  low  a  probability  of  deep  fades. 

Figure  7  compares  the  computations  with  the  two  model,  a*  X*  , 

0 . 1  and  s  =»  1 . 5  ( n  =  a  5  w  -  *  /  a  \  ,  Ry 

^  this  case,  the  computed 
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results  lie  midway  between  the  two  models,  with  Rice  underpre¬ 
dicting  and  Nakagami-m  overpredicting  the  probability  of  near- 
nominal  intensity  levels,  while  Nakagami-m  again  under  predicts 
and  Rice  now  somewhat  overpredicts  the  probability  of  deep  fades, 


simplest  parameterization  most  commonly  regarded  as  being  a 
nominal  overall  representation  of  the  ambient  ionosphere.  Of 
course,  it  must  be  recognized  that  no  single  PSD  parameteriza¬ 
tion  will  be  correct  in  general,  that  a  single  power  law  form 
may  often  be  too  simplistic,  and  also  that  conditions  in 
nuclear-disturbed  ionospheres  may  differ  appreciably  from  the 
ambient . 


Figure  8.  Comparison  of  computed  (C)  results  with  Rice  (R)  and 
Nakagami-m  (N)  models  at  s  =  2,  =  0.1,  s|  =  0.378, 


0.15b- 


0 . 0'5  — 


.01  .02  .  03  .  0*4  ,05 

Relative  Signal  Intensity  (I) 
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In  Figure  8  it  is  seen  that  at  s  =  2  both  Rice  and  Nakagami-m 
underpredict  the  probability  of  near-nominal  signal  intensity 
levels;  but  Rice  continues  to  overpredict,  and  Nakagami-m  to 
underpredict,  the  probability  of  deep  fades.  Comparison  of 
Figures  7  and  8  suggests  that  Nakagami-m  will  provide  a  good  fit 
to  the  computations  at  near-nominal  intensity  levels  for  a 
single  power  law  PSD  with  an  exponent  s  3  1.75  (n  3  3.75, 
v  £  11/8).  However,  the  deeper  fades  (although  relatively 
infrequent)  are  very  important  for  satellite  communications 
links;  and  even  at  s  =  2  Nakagami-m  continues  to  underpredict 
the  probability  of  fades  deeper  than  1?  dB  or  so  (i.e.,  at 


Lastly.  Figure  9  shows  the  corresponding  computed  results  ver¬ 
sus  model  predictions  at  =  0.1  for  s  =  2.5.  Here  it  is 

clear  that  both  models  tend  to  significantly  underpredict  the 
probability  of  near-nominal  signal  levels.  They  also  overpre¬ 
dict  the  probability  of  above-nominal  signal  enhancements;  a 
trend  towards  this  effect  can  also  be  perceived  in  the  previous 
figures.  Both  models  now  overpredict  the  occurrence  of  moder¬ 
ately  deep  fades,  Rice  much  more  so  than  Nakagami-m;  but  the 
Nakagami-m  model  continues  to  underpredict  the  likelihood  of 
fades  deeper  than  about  20  dB . 
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Figure  9. 


Comparison  of  computed  (C)  results  with  Rice  (R) 
Nakagami-ra  (N)  models  at  s  -  2.5,  X*  =»  0  1  s2 


and 
0.438 . 


At  other  values  of  X|y,  analogous  systematic  trends  are 

Observed  in  the  computed  signal  intensity  statistics  versus  the 
exponent  of  a  single  power  law  psn  anri  _ 

with  the  :  '  d  a°  ln  the  comparisons 

'  or  Nakagami-m  models.  The  details  of  these 
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trends  and  comparisons  will  differ  as  a  function  of  X*y  (see 
Appendixes  D  and  E).  The  general  inferences,  however,  remain 
as  indicated.  Rice  is  somewhat  understandably  f air-to-excel - 
lent  fit  at  s  <  3/2  and  generally  bounds  the  severity  of  the 
intensity  fluctuations  under  other  conditions.  Nakagami-m 
happens  to  be  a  useful  fit  near  s  3  2,  but  not  for  deep  fades. 
Lastly,  the  strong  and  systematic  trends  in  the  computed 
results  reveal  the  existence  of  signal  effects  that  cannot  be 
reproduced  in  general  by  any  simple  model  whicn  does  not 
account  for  the  details  of  the  PSD  parameterization  (and  also, 
as  demonstrated  shortly,  the  Fresnel  length). 

3-2  BER  CHARACTERISTICS. 

Analogous  trends  are  seen  in  the  computed  channel  symbol  BER- 
versus-SNR  characteristics.  Selected  examples  will  be  pre¬ 
sented  and  discussed  here,  using  DBPSK  as  a  point  of  reference. 
Numerous  additional  results,  including  other  modulations,  can 
be  found  in  Appendixes  E  and  F.  Similar  behavior  is  found  for 
all  modulations  considered. 

3.2.1  Variation  with  Spectral  Index. 

Figure  10  shows  the  computed  DBPSK  BER  characteristics  over  a 
range  of  values  and  for  two  extreme  values  of  the  PSD 

spectral  index  s.  The  X^y  values  are  listed.  Figure  10a  Is 
for  s=l  (n~2,  u=l),  and  Figure  10b  is  for  s  =  3  (n  *  5, 
u  =  2 ) .  The  overall  trend  is  apparent.  At  any  fixed  values 

■i4 

of  SNR  and  Xgy,  the  BERs  decrease  with  increasing  spectral  index. 
Similar  behavior  is  found  throughout  the  interval  1  <  s  <  3 
(see  Appendix  E). 
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Except  perhaps  at  very  low  levels  of  scattering,  these  effects 
are  manifested  even  more  strongly  when  the  BER  characteristics 


Figure  10.  Computed  DBPSK  channel  symbol  BERs  versus  SNR 

and  x]|y  for  single  power  law  PSDs  with:  (a)  s  *  1 
and  { b )  s  =  3 . 


are  parameterized  in  terms  of  rather  than  xjj[^.  This  is 
demonstrated  in  Figure  11.  On  the  left.  Figure  11a  (taken  from 
Fig.  3  of  Section  1)  has  fixed  s|  »  0.378  throughout;  the 
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corresponding  Xj^  values  decrease  with  increasing  spectral 
index  with  X^,  =  o.i  at  s  =  2.  On  the  right.  Figure  lib  has 
fixed  =  0.1  throughout.  The  corresponding  values  now 

increase  with  increasing  spectral  index;  but  despite  the 
increasing  value  of  the  BERs  at  fixed  x|y  and  SNR  actually 
decrease  with  increasing  s.  Neither  Nakagami-m  nor  Rice  would 
allow  such  an  effect  (i.e.,  a  decrease  in  BER  as  s|  is 
increased) .  It  is  apparent  that  the  reliability  in  detail  of 
simple  models  such  as  Rice  or  Nakagami-m,  which  are  specific¬ 
ally  parameterized  only  in  terms  of  s|,  can  be  very  dependent 
upon  the  actual  form  of  the  underlying  PSD. 


SNR  <CB)  SNR  <4Bi 

Figure  11.  Computed  DBPSK  channel  symbol  BERs  versus  spectral 
index,  either  at  fixed  »  0.378  (Fig.  11a)  or  at 
fixed  Xpy  =  0.1  (Fig.  lib). 
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3.2.2  Comparison  to  Rice  and  Nakagami-m. 


Additional  and  more  specific  comparisons  of  the  computed  BER 
results  with  Rice  and  Nakagami-m  statistics  are  presented  and 
discussed  here.  Only  four  examples  will  be  given;  many  more  can 
be  found  In  Appendix  F,  We  will  stress  intermediate  levels  of 
scattering  intensity,  specifically  =0.1,  in  order  to  most 

clearly  depict  the  differences  which  can  occur  between  the 
calculated  results  and  the  model  predictions.  At  much  smaller 
values  of  or  S£  all  results,  whether  from  the  computations 

or  from  the  models,  merge  together  towards  the  AWGN  limit;  for 
larger  or  1.0,  ©.11  results  merge  towards  the  SRF  limit. 

At  s  <  1,5  Rice  statistics  again  provide  a  fairly  good  fit  to 
the  computed  results  {Figure  12a).  At  s  =  2,  Nakagami-m  fits 
the  computed  results  fairly  well  at  SNR  £  15  dB  (Figure  I2b) 
but  underpredicts  the  BERs  it  higher  SNR  levels.  This  is  due 
to  the  tendency  of  Nakagami-m  to  consistently  underestimate  the 
frequency  of  deep  fades.  At  s  =  2.5,  this  particular  defici¬ 
ency  in  Nakagami-m  statistics  is  compensated  somewhat  by  an 
overprediction  of  the  more  shallow  fades  (see  earlier  Figure  9). 
Thus,  as  shown  in  Figure  12c,  Nakagami-ra  now  overpredicts  the 
computed  BERs  up  to  SNR  S  2i  dB ,  but  still  imderpredicts  the 
error  raters  at  higher  SNR  levels.  The  trend  continues,  as 
shown  1  r j  Figure  12d;  here,  at  s  =  3,  the  overprediction  by 
Nakagami-ra  of  the  shallower  fades  causes  a  corresponding  over¬ 
prediction  of  the  computed  BERs  up  to  SNR  3  30  dB,  At  higher 
SNR  levels,  Nakagami-m  would  again  underpredict  the  computed 
BER  levels. 
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It  is  clear  that  Rice  is  an  excellent  fit  at  s  =  1,  and  remains 
a  good  approximation  at  s  £  1,5.  For  higher  values  of  the  spec¬ 
tral  index  Rice  provides  an  upper  bound  to  the  computed  BERs . 

For  reasons  discussed  above,  Nakagami-m  never  reproduces  the 
computed  BER-versus-SNR  characteristic  in  detail,  but  happens 
nonetheless  to  provide  a  good  fit  at  s  *  2  and  SNR  £  15  dB, 
or  an  approximate  fit  at  s  =  2.5  and  SNR  £  21  dB . 

3.3  THE  CASE  sj  >  1.0. 

There  is  one  special  case  which  merits  at  least  brief  discus¬ 
sion  in  the  context  of  single  power  law  PSD  parameterizations. 
This  is  the  case  s|  >_  1.0,  which  is  not  allowed  by  either  Rice 
or  Rayleigh  statistics,  nor  {It  would  seem)  purposefully  accom¬ 
modated  In  the  design  of  the  Nakagaml-m  model. 

With  signals  of  fluctuating  intensity,  >  l.o  merely  implies 
that  the  standard  deviation  of  the  signal  intensity  is  greater 
than  its  mean  value.  Such  a  thing  is  clearly  possible,  and  our 
specific  computations  also  show  that  conditions  of  SJ  >  1.0  can 
easily  occur  with  PSD  parameterizations  encompassing  those  of 
current  interest.  There  might  be  justifiable  concern  that  s|  > 
1.0  could  lead  to  BER  levels  which  might  greatly  exceed  the 
normally  accepted  "SRF  limit"  for  slow  fading  effects.  As  a 
practical  matter,  however,  our  computations  suggest  that  the 
resulting  slow  fading  BER  characteristics — over  the  range  of 
PSD  parameterizations  investigated--aiways  J  te  either  below  or 
within  a  few  dB  or  less  above  the  SRF  limit. 

To  depict  the  general  behavior  encountered,  we  will  first  present 
two  examples  in  Figure  13,  below;  many  more  can  be  found  in 
Appendix  C.  As  already  shown  in  Table  1  and  Figure  4.  a  signifi¬ 
cant  tendency  towards  conditions  with  s|  >  1,0  clearly  exists  at 
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s  £  2,  but  it  is  most  pronounced  at  the  higher  values  of  spectral 
index.  Therefore,  we  will  concentrate  here  on  the  casts  3  =  3. 

Figure  13  depicts  most  of  the  features  of  the  anomalous  s|  > 

1.0  case  as  encountered  in  our  computed  results.  The  probabil¬ 
ity  of  very  large  signal  intensity  enhancements  tends  to  sig¬ 
nificantly  exceed  that  predicted  by  Rayleigh  statistics.  Yet, 
for  more  moderate  levels  of  s  gnal  enhancement,  the  computed 
intensity  distributions  lie  systematically  below  the  SRF  limit. 
Then,  at  slightly  sub-nominal  signal  intensity  levels  (i.e., 
very  weak  fades) ,  the  computed  intensity  distributions  again 
exceed  the  3  1.0  predictions  of  Rayleigh  and  other  models. 

The  phenomena  described  above  are  observed  consistently  in  all 
our  computed  results  for  s|  >  1.0.  At  lower  levels  of  signal 
intensity,  however,  a  transition  is  observed.  As  shown  in 
Figure  13a,  at  moderate  levels  of  X^y  the  computed  probability 
of  deep  fades  is  less  than  that  predicted  by  Rayleigh  statistics. 
At  higher  levels  of  X^y,  the  computed  results  gradually  merge 
toward  the  Rayleigh  case.  Throughout  this  interval,  however, 
the  computed  channel  symbol  BKRs  continue  to  lie  within  a  few 
dB  or  less  of  (but  sometimes  above)  the  "SRF  limit." 

Typical  BER  behavior  at  >  1.0  is  illustrated  in  Figure  14, 
which  contains  the  same  two  cases  from  Figvre  13,  plus  a  third 
example  with  X^y  »  0.25  and  =*  1.28.  For  X^y  =  0.25  to  0.40, 
and  =  1.28  to  1.65,  the  BER  characteristic  still  lies  appre¬ 
ciably  below  the  SRF  limit,  just  as  it  does  for  other  PSDs  with 
s  >  1.5  spectral  indexes  at  these  moderate  or  lower  levels  of 
Xpy.  For  X^y  »  7.0,  the  value  of  is  1.99:  this  is  essenti¬ 
ally  twice  the  value  { *  1.0)  associated  with  Rayleigh  sta¬ 
tistics,  and  it  is  also  nearly  the  largest  value  observed  in 
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Relative  Signal  Intensity  (I|  Relative  Signal  Intensity  (I) 


Figure  13.  Comparison  of  computed  signal  Intensity  distribu¬ 
tions  for  Sj  >  1 . 0  at  s  *  3  with  Rayleigh  statis¬ 
tics:  (a)  X*y  -  0.4,  (b)  xfy  =.  7.0. 

all  of  our  computed  results  (see  Table  1  or  Figure  4).  None¬ 
theless,  the  computed  BER  characteristic  never  exceeds  the  SRF 
limit  by  more  than  about  2  dB;  and  it  will  then  move  closer  to 
SRF  as  X|y  increases.  Thus,  despite  the  unusual  and  highly  non- 
Rayleigh  oehavlor  of  the  computed  Intensity  distributions  at  > 
1.0,  the  corresponding  BER  characteristics  never  substantially 
exceed  these  which  would  be  predicted  by  Rayleigh  statistics. 
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SECTION  4 

TWO  POWER  LAW  RESULTS 


Our  computed  results  for  first-order  signal  intensity  statis¬ 
tics  with  a  typical  and  realistically  parameterized  "two  power 
law"  PSD  were  found  to  be  entirely  analogous  to,  and  qualita¬ 
tively  explainable  by,  the  results  from  single  power  law  PSDs 
of  varying  spectral  index,  as  already  presented  in  Section  3. 
Therefore,  a  smaller  number  of  specific  examples  for  two  power 
lav/  PSDs  will  be  presented  and  discussed  here.  Additional 
examples  can  be  found  in  the  Appendixes. 

The  representative  form  of  the  two  power  law  PSD  that  we  employ 
is  given  in  Equation  6  of  Section  2,  with  the  constants  deter¬ 
mined  by  Equation  A-6  of  Appendix  A.  The  selected  values  for 
v2  and  kB  are: 

Xj  =  5/4 ,  (Sj  =  1 . 5) , 

v2  =  2,  (S2  =  3 ) , 

kg  =  2n/ 700  m” 1 . 

These  are  representative  of  typical  values  determined  from  fit¬ 
ting  a  more  complex,  two  power  law,  PSD  to  ambient  environment 
data  (Ref.  37).  This  PSD  form  is  shown  in  Figure  15.  The 
Fresnel  lengths  (Equation  13)  which  we  will  consider  in  great¬ 
est  detail  are  (in  meters) 

*F(«)  =  1386,  693,  346,  173,  «  =  1-4, 

and  their  locations  in  k-space,  k# ,  are  also  shown  in  Figure  15 
where  kR  =  2n/Jp(«). 

We  might  anticipate  on  the  basis  of  Figure  15  and  earlier  discus 
sions  that  the  kj  results  would  be  mainly  character 1st ic  of  the 
Sj  slope  and  the  k4  results  would  be  more  nearly  characteristic 
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of  the  s2  slope.  The  k2  and  k3  results  should  be  Intermediate, 
although  there  will  be  a  tendency  toward  the  e.^  slope;  this  is 
expected  since  the  results  really  depend  on  a  region  of  k-space 
that  includes  small  v.  Lues  of  k  up  through  and  somewhat  beyond 
the  vicinity  of  ka.  Consequently,  even  the  k3  and  k4  results 
will  still  be  influenced  by  the  slope  region.  In  Section  1, 
we  very  loosly  summarized  this  complex  behavior  by  indicating 
that  the  results  are  roughly  equivalent  to  those  of  a  single 
power  law  with  an  effective  slope  determined  by  the  behavior  at 
wave-vectors  "generally  less  than  and  in  the  vicinity  of"  kK „ 


Figure  15.  Two  power  law  power  spectral  density;  dso  indicated 
are  the  locations  of  ka  («  »  1-5). 
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In  addition  to  the  detailed  results  discussed  below  for  *  *  1-4, 
we  will  also  present  a  few  results  for  «  -  5  with 

Ap  (  5  )  =  ?  7  m , 

that  is,  a  value  of  kw  twice  as  large  as  k^ ,  and  thus  further 
removed  for  the  spectral  break.  As  anticipated,  this  Fresnel 
length  leads  to  behavior  ever:  more  characteristic  of  Sg,  but 
the  influence  of  is  still  evident. 

In  this  section,  we  will  discuss  the  trends  as  the  Fresnel 
length  is  changed.  These  trends  will  also  be  illustrated  by  a 
few  examples.  Additional  results  can  be  found  in  the  Appen¬ 
dixes.  Appendix  G  presents  tha  calculated  intensity  prob¬ 
ability  distributions.  The  small  intensity  behavior  is  given 
in  Appendix  H.  Since  the  main  conclusions  regarding  small 
intensity  behavior  are  essentially  the  same  as  for  the  corres¬ 
ponding  single  power  law  PSD3 ,  we  will  not  discuss  them  further 
in  this  section.  The  channel  symbol  BER  curves  are  found  in 
Appendix  I,  and  the  comparison  with  Rice  and  Nakagami-m  for 
DBPSK  is  in  Appendix  J. 

4.1  SIGNAL  CHARACTERISTICS. 

As  with  the  single  power  law  PSD,  one  of  the  parameters  that 
characterize  the  probability  distributions  is  the  scintillation 
index  sf.  Equation  13.  Table  2  presents  a  complete  list  of  the 
calculated  values  of  S^  for  the  selected  values  of  and  ka. 

Including  «  *  5.  This  same  information  is  presented  graphic¬ 
ally  in  Figure  16.  For  small  values  of  X^y(<  0.1),  all  s| 
versus  X^y  results  continue  to  approximately  satisfy  the  linear 
relationship,  S*  2  4Xj|y .  As  X^y  increases,  and  for  the  Fresnel 
lengths  larger  than  or  on  the  order  of  the  freezing  length  («  * 
1,2),  s|  overshoots  and  then  returns  towards  unity.  For  she 
shorter  Fresnel  lengths,  considerable  transitory  overshoot  of 
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is  possible.  As  we  saw  in  the  previous  section,  such  over¬ 
shoot  of  unity  for  s|  is  characteristic  of  steep  (s  >  2)  PSDs . 
However,  as  also  anticipated,  the  overshoot  is  somewhat  mod¬ 
erated  here  by  the  influence  of  the  s  *  1.5  PSD  behavior  at 
small  values  of  k. 


The  signal  intensity  distributions  also  show  the  trend  from 
approximately  2  »  1 . 5  behavior  toward  approximately  s  »  3 
behavior  as  k^  increases.  For  5^  <  1 ,  the  distributions  tend 
to  be  fairly  close  to  Nakagami-m  at  moderate  intensity  levels 
and  for  the  longer  Fresnel  lengths  considered  («  *  1,  2). 
However,  Nakagami-m  continues  to  underestimate  the  probabii  ity 
of  fades  deeper  than  10-to-20  dB  (see  deep  fade  results  in 
Appendix  H) ;  and  Rice  continues  to  overestimate  these  fades. 


Table  2 . 


vs 


S 


2 

4 


and  kK. 


CASE 

X* 

Rv 

«  =  1 

«  *  2 

«  «  3 

*  -  4 

«  =  5 

a 

0 .01 

0.040 

0.039 

0.040 

0.041 

0.041 

b 

0.025 

0.097 

0.097 

0.103 

0.107 

0.109 

c 

0 . 05 

0.186 

0 . 188 

0.208 

0.227 

0.238 

d 

0.1 

0.347 

0.350 

0.409 

0.479 

0.519 

e 

0.25 

0.672 

0 . 684 

0.827 

1 . 02 

1  .  17 

f 

0.4 

0.837 

0.854 

1.02 

1 . 25 

1 . 43 

g 

0.7 

C.965 

0.S98 

1  .  17 

1 . 41 

1.65 

h 

1 . 5 

1.03 

1 . 09 

1 . 23 

1.45 

1.70 

i 

3.0 

1 .05 

1.11 

1 . 23 

1 . 41 

1.63 

J 

7.0 

1.05 

1  .  10 

1 . 20 

1 . 32 

1.51 

k 

10.0 

1 . 04 

i  .  10 

.1 . 18 

1 . 28 

1 .46 

1 

20.0 

1.03 

1.08 

1  .  14 

1.23 

1.37 
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As  the  Fresnel  .length  becomes  shorter  than  the  freezing  length 
{*  =:  3,4),  even  the  limited  agreement  with  Nakagomi-m  at  near- 
nominal  signal  levels  disappears.  As  examples,  Figure  17  shows 
the  calculated  distributions  for  X^y.  =  0.05,  «  *  1  and  4. 


Figure  17.  Calculated  signal  intensity  distribution  for  two 

power  law  striation  PSD  of  Figure  15,  at  X*  =0.05 

«y 

for  Ap  *  13P6  m  and  4p  =  173  ra,  as  compared  to 
Rice  and  Nakagami-m. 
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The  intensity  distributions  for  S4  >  1  also  show  the  antici¬ 
pated  trends  based  upon  the  single  power  law  PSD  results.  As 
the  Fresnel  length  is  decreased  the  expected  progression  from 
Rayleigh  to  distributions  that  are  larger  than  Rayleigh  for 
both  small  and  large  values  of  the  intensity  is  observed. 
Detailed  results  are  given  in  Appendix  G. 

4.2  CHANNEL  SYMBOL  BIT  ERROR.  RATES. 

As  with  the  single  power  law  results,  0.01  $  X^y  5t  0.4  essen¬ 
tially  spans  the  gap  between  the  AWGN  and  SRF  limits  for  the 
BER  curves.  Likewise,  the  results  for  X^,  =  0.1  are  of  partic¬ 
ular  interest  and  are  shown  in  Figure  18  (including  «  =  5)  for 
DBPSK,  This  figure  again  illustrates  that  X^y  “0.1  is  a 
reasonable  compromise  between  the  AWGN  and  SRF  limits.  As 
would  be  anticipated  from  the  single  power  law  results,  the 
BER  improves  as  the  Fresnel  length  decreases  for  this  fixed 

value  of  .  Figure  18  can  be  compared  with  Figure  lib  of 
Ry 

Section  3,  Recall  (Fig.  11a)  that  a  larger  spread  in  BER 
characteristics  will  be  seen  at  fixed  S4  than  at  fixed 

The  comparisons  of  the  BER  results  with  Rice  and  Nakagami-m  are 
also  as  expected.  For  *  =  1  the  calculated  results  agree 
reasonably  well  with  Rice.  At  e  =  2 ,  they  may  be  compared 
with  Nakagami-m  except  at  large  SNR  where  Nakagami-m  again 
underestimates  the  BER.  For  «  =  3 ,  4,  the  BER  curves  behave 
quite  differently  from  either  Rice  or  Nakagami-m.  Rice  always 
overestimates  the  BER,  and  typically  there  is  a  range  of  lower 
SNR  levels  for  which  Nakagami-m  also  overestimates  the  BER.  As 
illustrations ,  the  comparisons  at  =  0.1  for  «  =  1,  5  are 

shown  in  Figure  IS  For  larger  values  of  xj|y  (Xpy  )  0.4),  the 
BERs  are  usually  quite  close  to  SRF.  Small  deviations  are 
observed,  which  are  very  similar  to  those  already  discussed  in 
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Section  3.  As  expected,  the  deviations  are  greatest  for  »  kg. 
Figure  19  may  be  compared  with  Figure  12  of  Section  2.  Thus, 
Figure  19a  at  «  *  1  shows  behavior  nearly  identical  to  that  of 
Figure  12a  at  s  »  1.5;  and  Figure  19b  at  «  =  5  is  moat  similar 
to  Figure  12d  at  s  «  3. 


SNR  (dB) 

Figure  18.  BER  vs.  SNR  characteristics  for  DBPSK  for 

various  Fresnel  lengths  Jp(«)  at  »  0.1. 


1 


SNR  (dB)  SNR  fdB) 


Piguro  19.  Calculated  channel  symbol  bit  error  rate  vs.  SNR 
at  =  0.1  and  i p  =»  1386  m  and  ip  *  86  m,  as 

compared  to  Rice  and  Nakagami-m. 
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SECTION  5 
DISCUSSION 


The  final  discussion  will  consist  of  two  parts.  Firstly,  we 
will  supplement  the  preceding  computational  results  with  a 
review  of  past  analyses  and  findings  from  the  pertinent  data  on 
satellite  link  intensity  distributions  in  the  ambient  environ¬ 
ment.  We  will  then  conclude  with  a  brief  summary  of  what  has 
been  learned  and  its  implications  for  further  research  and 
systems  applications. 

5.1  REVIEW  OF  AMBIENT  ENVIRONMENT  DATA. 

Whitney,  et  al .  (Ref,  32)  published  one  of  the  earliest 
experimental  studies  of  satellite  link  first-order  intensity 
statistics.  They  compared  their  data  only  with  Nakagami-ra. 

The  data  were  grouped  together  into  five  categories  according 
to  the  value  of 


SI  S  100 


Xmax  Jmin 
^■max  +  Imin 


For  each  group,  an  overall  median  intensity  distribution  was 

determined  from  the  data  and  then  compared  to  Nakagaml-m.  From 

the  comparison,  a  best-fit  value  for  s|  was  determined.  Thus, 

2 

S4  was  actually  treated  as  an  adjustable  parameter.  With  this 
degree  of  adjustment  at  their  disposal,  Whitney,  et  al  were 
able  to  find  five  different  s|  values  which  would  allow 
Nakagaml-m  to  fit  the  overall  median  intensity  distributions  of 
each  of  their  five  different  SI  groups. 


The  fit  to  individual  data  sets  from  each  SI  group  was  poorer. 
The  fit  was  also  found  to  deteriorate  noticeably  at  fade  depths 
on  the  order  of  ~10  dB  or  greater;  and  very  few  data  were 
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available  at  fade  depths  beyond  ~10-to-lfi  dB  or  so.  Nonethe¬ 
less,  Whitney  et  al  recommended  that  Nakagami-m  could  be  used 
reliably  to  extrapolate  from  the  data  to  greater  depths  of 
fade.  Based  upon  our  present  computed  results,  it  seems  that 
this  conclusion  was  not  correct;  even  under  conditions  where 
Nakagami-ra  gives  a  good  fit  for  shallow  fades,  it  does  not 
accurately  match  deep  fade  statistics. 

In  a  subsequent  study,  Rino  and  Fremouw  (Ref.  33)  tested 
several  hypothesized  simple  models  against  data.  Considered 
were:  Nakagami-m,  Rice,  Gaussian,  and  log -normal  statistics. 

They  concluded  that  equally  good  fits  could  be  obtained  from 
either  the  Gaussian  or  log-normal  distribution,  but  that  either 
Rice  or  Nakagami-m  gave  only  a  poor  approximation  to  the  data. 

The  data  used  in  this  study  were  also  of  limited  dynamic  range, 
and  therefore  did  not  accurately  sample  deep  fades.  Nonetheless, 
the  conclusions  reached  by  Rino  and  Fremouw  were  clearly  very 
different  from  those  of  Whitney,  et  al . 

In  yet  a  subsequent  study,  Rino,  et  al  (Ref.  30)  compared 
additional  data  sets  with  Rice  statistics  and  with  generalized 
Gaussian  and  log-normal  models.  According  to  a  later  report  by 
Fremouw,  et  al  (Ref.  31 — see  below),  these  authors  did  not 
separately  consider  Nakagaml-m  "because  they  viewed  it  as 
virtually  Identical  to  the  Rice  distribution,"  although  It 
seems  difficult  at  present  to  Imagine  that  such  a  misperception 
could  really  have  occurred.  In  any  event,  Rino,  et  al  con¬ 
cluded  that  Rice  gave  a  rather  poor  fit  to  their  measured  data, 
but  that  a  generalized  Gaussian  (with  two  adjustable  parame¬ 
ters)  gave  a  somewhat  better  fit  than  the  corresponding  log¬ 
normal  form  (with  no  free  parameters). 
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Lastly,.  Fremouw,  et:  al  (Ref.  31)  performed  what  is  apparently 
the  most  extensive  model-data  comparison  of  note,  and  reached 
yet  a  different  set  of  conclusions.  In  this  study,  Nakagami-m, 
log-normal,  and  generalized  Gaussian  models  were  tested  against 
the  data,  as  well  as  a  two-component  model  in  which  the  signal 
is  considered  to  be  the  product  of  a  focused  component  with 
log-normal  statistics  and  a  scattered  component  with  generalized 
Gaussian  statistics.  Except  for  Nakagami-m,  the  models  were 
also  tested  for  their  ability  to  fit  signal  phase  statistics, 
as  well  as  intensity  statistics. 

The  basic  conclusion  of  Fremouw,  et  al  was  that  Nakagami-m  gave 
the  better  fit,  in  agreement  with  the  original  finding  of 
Whitney,  et  al . ,  but  in  contradiction  to  the  more  recent  con¬ 
clusions  of  Rino  and  Fremouw  and  (  inl'erentially  at  least)  of 
Rino,  et  al ,  More  importantly,  however,  the  fit  of  Nakagami-m 
to  the  data  (which  still  did  not  accurately  sample  deep  fades) 
was  simply  not  very  good.  That  is,  Nakagaml-m  gave  the  best 
fit  in  only  32  out  of  83  cases,  while  generalized  Gaussian  did 
nearly  as  well  as  best  fit;  and  one  or  the  other  of  the  two 
remaining  models  also  gave  a  best  fit  more  that  27  percent  of 
the  time. 

In  hindsight,  the  proper  conclusion  to  be  drawn  from  all  this 
is  simply  that  none  of  these  models  are  reliable  in  general, 
even  for  only  the  ambient  environment,  and  even  when  deep  fades 
are  poorly  sampled.  The  cause  is  evident  in  our  computed 
results.  The  actual  signal  intensity  distributions  will  have  a 
strong  and  systematic  dependence  on  the  power  spectrum  of  the 
scattering  ionospheric  structure  as  well  as  the  Fresenl  length 
in  general;  this  structure  (and  the  wavelength-dependent  Fresnel 
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.length)  is  at  least  somewhat  variable,  even  for  the  ambient 
environment;  and  none  of  the  simple  models  put  forth  to  date 
have  beer,  designed  to  accommodate  these  effects. 

Moreover,  the  fact  that  Nakagami-m  is  not  inf requently  a  good 
fit  to  the  ambient  environment  data  when  deep  fades  are  not 
included  is  entirely  consistent  with  our  computed  results, 
since  we  also  find  Nakagami-ra  to  be  a  useful  heuristic  fit  at 
moderate  intensity  levels  for  the  PSD  conditions  (i.e.,  «»  k~3 
form  for  k  £  j»p)  thought  to  most  nominally  apply  for  the 
ambient  environment. 

5.2  CONCLUSIONS  AND  OUTLOOK. 

The  many  computed  results  presented  In  this  report  and  Its 
Appendixes  clearly  reveal  that  signal  intensity  statistics  in 
non-Rayleigh  fading,  as  well  as  the  basic  regime  of  validity  of 
the  Rayleigh  fading  approximation,  are  strongly  dependent  upon 
the  parametric  representation  of  the  scattering  medium  and 
the  value  of  the  Fresnel  length.  Although  the  specific  depend¬ 
encies  that  we  have  demonstrated  have  been  based  on  the  thin 
phase  screen  approximation  and  its  details,  it  is  clear  that 
they  are  of  mere  fundamental  origin.  Since  none  of  the  exist¬ 
ing,  simple  models  for  non-Rayleigh  fading  have  been  designed 
with  these  effects  in  mind,  it  is  not  surprising  that  such 
models  are  not  reliable  in  general. 

The  approximate  agreement  of  Nakagami-m  with  either  data  or 
computations  within  certain  ranges  of  PSD  and  Fresnel  length 
parameterization,  and  also  within  certain  regimes  of  signal 
Intensity  fluctuation,  appears  to  be  at  least  somewhat  a  coin¬ 
cidence,  and  of  nonphysical  origin.  Nonetheless,  Nakagami-m 
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does  happen  to  provide  a  useful  fit  for  nominal  ambient 
environment  PSDs  and  for  fade  depths  or  average  SNR  levels 
which  are  not  too  large. 

Rice  intensity  statistics  are  thought  to  be  more  physically 
based  and  they  provide  a  very  good  fit  oc  all  levels  of  signal 
intensity  or  SNR  for  conditions  which  can  be  represented  by  a 
PSD  with  spectral  index  s  <  1.5.  However,  these  physically 
definable  conditions  are  apparently  more  stressing  than  com¬ 
monly  encountered  in  the  ambient  ionosphere;  and  it  appears, 
with  some  uncertainty,  that  they  may  also  represent  an  extreme 
case  for  the  nuclear-perturbed  environment. 

Thus,  in  general,  it  appears  that  Nakagami-vn  can  continue  to  be 
used  as  a  convenient  but  largely  empirical  curve  fit  for  ambient 
environment  conditions  at  link  SNRs  below  15  dB  or  so,  while 
Rice  can  be  used  as  a  more  fundamental  approximate  upper  bound 
on  the  severity  of  the  signal  fading  under  all  ambient  or 
nuclear-perturbed  conditions  of  probable  interest.  The  useful¬ 
ness  of  an  empirical  fit  of  limited  validity  (i.e.,  Nakagami-m) 
versus  a  potentially  more  rigorous  upper  bound  (i.e..  Rice) 
will  depend  upon  the  priorities  of  the  application. 

For  more  precise  work,  the  alternative  in  principle  is  to  use 
detailed  computations  such  as  those  employed  in  the  present 
study,  or  even  as  generalized  to  bypass  the  thin  phase  screen 
approximation.  However,  the  credibility  of  this  more  laborious 
alternative  depends  entirely  upon  the  ability  to  specify  the 
phase  screen  PSD  or  its  underlying  refractive  index  power 
spectrum  in  an  appropriate  level  of  detail.  The  reliability 
with  which  this  environment  specification  can  now  or  eventually 
be  provided  is  largely  a  matter  of  judgement  and  current 
research.  In  applications  which  require,  and  can  afford  the 
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cost  of  insensitivity  to  such  uncertainties.  Rice  statistics 
will  continue  to  provide  a  useful  "worst  case"  specification 
for  the  signal  intensity  statistics  under  all  levels  of 
scattering  intensity,  including  the  Rayleigh  fading  limit. 

Our  findings  from  this  investigation  indicate  two  priorities 
for  further  research.  On  the  one  hand,  research  aimed  at 
better  and  more  reliable  characterization  of  the  scattering 
medium  must  clearly  seek  to  reduce  the  environment  of  uncer¬ 
tainties  to  the  point  that  they  no  longer  dominate  the  uncer¬ 
tainty  in  intensity  statistics  for  non-Rayleigh  fading,  or 
the  uncertainty  in  the  value  of  the  election  density  or 
refractive  index  variance  needed  to  drive  the  weak- to-strong 
scattering  transition  (and,  thus,  the  regime  of  validity  of 
Rayleigh  statistics). 

On  the  other  hand,  the  fact  that  the  Rice  distribution  seems  to 
offer  a  useful  "worst  case"  specification  for  first-order 
intensity  statistics  suggests  the  possibility  that  a  generali¬ 
zation  of  Rice  intensity  statistics  to  provide  a  reasonably 
worst-case  specification  for  both  amplitude  and  phase,  as  well 
as  for  second-  and  higher-order  signal  statistics,  may  also  be 
achievable.  This  is  an  attractive  concept,  since:  (1)  Rice 
statistics  merge  to  Rayleigh  statistics  in  the  strong  scattering 
limit;  (2)  both  Rice-  and  Rayleigh- type  signal  realizations  are 
easily  generated  from  sampling  Gaussian  statistics;  and  (3)  the 
alternative  of  using  Presnel-Kirchhof f  calculations  is  both 
difficult  and  subject  to  the  reliability  with  which  the  under¬ 
lying  refractive  index  spectrum  can  be  defined. 

One  problem  to  be  met  in  pursuit  of  this  last  objective  is  the 
fact  that  Rice  statistics  will  not  be  "accurate"  in  general. 
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Instead,  they  will  presumably  be  an  approximate  worst  case;  but 
this  means  that  the  proper  terms  of  reference  for  judging  the 
utility  of  various  alternate  Rice  model  generalizations  remains 
somewhat  to  be  determined.  A  further  problem  to  be  met  is  the 
proper  separation  of  diffractive  and  refractive  signal  effects 
in  a  Rice-type  model.  Although  this  issue  most  likely  will  not 
affect  link  performance  assessments  under  those  weak,  slow 
fading  conditions  for  which  non-Rayleigh  signal  statistics  will 
pertain,  the  proper  specification  of  nondif f ractive  phase 
effects  in  general  is  an  open  issue  which  would  inevitably 
merit  further  consideration  during  the  development  of  a  "worst 
case"  signal  structure  specification  encompassing  both  Rayleigh 
and  non-Rayleigh  fading. 

In  the  final  analysis,  the  need  and  utility  of  such  further 
theoretical  work  towards  signal  structure  specifications  for 
non-Rayleigh  fading  will  depend  upon:  (1)  the  reliability  with 
which  ionospheric  structure  power  spectra  can  be  specified,  for 
both  ambient  and  nuclear-perturbed  environments  ( including 
multi-burst  nuclear  environments);  (2)  the  assessed  capability 
and  readiness  of  communications  system  design,  assessment,  and 
test  personnel  to  use  detailed  Fresnel-Kirchhof f  calculations 
versus  worst-case  signal  structure  specifications;  and  (3)  the 
probable  sensitivity  of  DoD  satellite  communications  systems  and 
networks  either  to  link  conditions,  including  non-Rayleigh 
fading,  or  the  reliable  definition  of  the  extent  of  the 
Rayleigh  fading  regime . 
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APPENDIX  A 

ADDITIONAL  MATHEMATICAL  DETAILS 


Section  2  presented  the  physical  and  mathematical  foundations 
of  tne  calculations  discussed  in  this  report.  We  here  extend 
and  comment  on  some  of  the  material  given  there. 

A . 1  SCINTILLATION  INDEX  AND  RYTOV  PARAMETER. 

Q  A 

The  relationship  between  and  Xp  for  weak  scattering 

(S|  ~  4  Xpy)  can  be  derived  as  follows.  From  the  expression 

for  sj.  Equation  13,  and  scattered  field.  Equation  9,  we  have 

s4  +  1  =  7777^2  Jdx1dx2dx3dx4  e*P  [17P’(X1'5<2+X3"X4)] 

X<eit'*l'*2++3"*A'>.  (  A- 1  ) 

In  the  limit  of  weak  scattering  and  employing  various  symmetry 
properties  of  the  under  the  sign  of  integration,  we  find, 

see  Equation  1, 


<ei(^l  ^2 +  ^3  -*  1-  I  ^4B(  0)  -8B(x1-x2 


+2B ( X ^ -x^ ) +2B ( Xj-x^ ) j .  ( A- 2 ) 


Next,  noting  that 


Jdx  exp  [Vi  x2J  =  [\Z*]**  e±ln/4,  (A-3) 


6  1 


we  easily  see  that  the  terras  independent  of  8  cancel  and 


sj  tr  4 

1 

XZ* 

Jdx1dx2  B(x1-X2)  |  exp  j 

[i  A 

•  XZ 

(xf-x^)] 

+  i 

pyn 

[*  IP(X^  +  X2>^  *  *  *XP 

fi  JL. 

/ v2+X« ) 1  1 

4 

CA(v 

1  1  * 
1  xz 

ix1+x2;j  > 

-  2B  (  0  )  . 

a 

Finally,  introducing  the  PSD,  Equation  2,  and  performing  the 


coordinate  integrations,  we 

have 

s?  -  4  1  f;  [i-  j  «*p  1 

[-i  Hi  k2) 

1  2n  f 

1 

-  -  exp 

(*  £  *2> 

1 

2 

5  4  I  r  +(k)  sin2  k2 

0 


=  4 


(A-5) 


A. 2  TWO  POWER  LAW  COEFFICIENTS. 

The  form  of  the  two  power  law  spectrum  that  we  have  used  is 
given  in  Equation  6.  The  coefficients  are  determined  by  con¬ 
tinuity, 

S2v2  Cl  =  c2(l  +  «2)*'l, 

and  normalization 

n  VS  T(  1  S 

Cj  2  nifj)  2  v»2  - 1  (1+S2)V1 


( A-6a) 


J 


dt 


i  ( i +t 2 ) 


( A-6b) 
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Here  we  have  defined  $  as 


and  is  typically  much  larger  than  unity, 

A.  3  RYTOV  PARAMETER  EVALUATION „ 


The  Rytov  parameter  was  defined  In  Equation  14. 


of  i , 


In  the  lire! t 

(A-  8) 


being  small,  which  is  the  typical  case,  we  have  for  the  single 
power  law  spectrum  (1/2  <  .  <  5/2) 


<n/2)* 


_ Sly) _ 

r (v-iTTTv+fc) 


(2 


n  n 

sm  -  cos 


(A '-9  ) 


Of  course,  the  v  =  1  spectrum  that  we  used  was  not  exactly  a 
strict  power  law,  Equation  5,  but  Equation  A-9  is  still  reason¬ 
ably  accurate.  However,  in  the  numerical  calculations ,  we 
actually  integrated  Equation  14  for  the  v  =  1  case.  In  addi¬ 
tion,  this  integration  was,  ot  necessity,  required  as  well  for 
the  two  power  law  PSD. 

A. 4  RELATIONSHIP  TO  ELECTRON  DENSITY  VARIANCE. 


It  was  noted  in  Section  2  that  there  is  a  connection  between 

at,  the  variance  of  the  LOS  phase  change,  and  the  variance  of 
■  / 
the  electron  density  fluctuations,  a£e .  That  connection,  for 

a  single  power  law  PSD  with  v  >  1,  is  given  in  Equation  12. 


The  corresponding  connection  for  v  =  1  is 


2  2 


r;  L  Lq  orNe  3— ——7  exp  i-ii/LQ) 


K0{JVLo) 


-  nX2  r  2  L  Lc  *fie/  Jn(  LQ/Jti ) 


A- 10 


and  for  the  two  power  law  spectrum,  Equation  6,  is 


«r|  =  2n  X2  r2  L  LQ  °rNe/c3' 


;  A-lla) 


3  *  5  *  |  jpi  *2<1-*'3,c2-  (A-Xlb) 


A. 5  BER  FOR  AWGN . 


In  the  text,  DBPSK  was  the  modem  selected  to  illustrate  the 
effects  of  various  levels  of  fading.  In  all,  we  have  con¬ 
sidered  the  behavior  of  six  modems.  The  BER 1 s  for  AWGN  are 
(Ref.  39) 


CPSK : 


aPSK: 


DBPSK: 


Pe  *  \  erfc  'H 

Pe  =  erfc  /y  (i-  ~  erfc  /r) , 


(A-12a) 

( A-12b) 


(A-120 


BFSK:  Pfe  =  i  e  1/2  , 

( A--12d ) 

QFSK:  Pe  =  1  -  [ 1  “  \  « 

"  T  ]  3  /  2  ^ 

( A-l 2e ) 

8  --ARYFS’C  :  P_  =  l-  [l--£e 

-3r/2j7/3( 

( A- 1 2  f ) 

where  g  is  the  bit  energy- to-noise  ratio  ,'SNR). 

a  fading  environment  is  given  in  Equation  20 

The  BER  in 

00 

<Pe(r)>  *  Jdi  P(i) 

Pe  (YD- 

(A-13) 

o 


It  can  be  seen  that  Pe  drops  off  very  rapidly  so  that  for  large 
values  of  y  in  Equation  A-13,  only  small  values  of  I  contribute. 
Consequently,  the  small  intensity  (deep  fades)  behavior  of  the 
probability  distributions  control  the  BER ' s  for  large  SNR.  If 
the  lading  statistics  is  Rayleigh,  the  corresponding  BER  is  the 
slow  Rayleigh  fading  limit  (SRF). 
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APPENDIX  B 

NUMERICAL  CONSIDERATIONS 


B.l  PHASE  SCREEN. 

The  realization  of  the  phase  screen  j.s  given  as  a  FFT  in 
Equation  8.  It  is  clear  from  its  creation  that  if  we  extend 
the  phase  screen  beyond  the  range  0(4^  N-i,  it  is  periodic, 

f*+N  ”  »'  (B-l) 

and  continuous.  As  an  illustration  of  this  last  point, 

Figure  20  shows  a  section  of  a  representative  normalized 
phase  screen,  generated  using  Equation  8  with  v  **  3/2 

and  N  =  16384.  We  here  plot  the  region 

A  -  15500  -  17500  (Mod  N)  (B-2) 

to  emphasize  that  there  is  no  discontinuity  at  J  *  H.  Such 
a  discontinuity  would  have  lead  to  edge  diffraction  when  the 
scattered  field  was  calculated, 

B.2  SAMPLED  FRESNEL-KIRCHHOFF . 

The  sampled  version  of  the  Fresnel-Kirchhof f  integral  was  given 
in  Equation  11.  This  result  depended  on  the  fact  that  the  main 
contribution  to  h{x),  Equation  9,  comes  from  the  region  of  x' 
near  x  and  we  have  sampled  x  and  x'  as 

x  »  max , 

X'  *  iAX. 
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Figure  20.  Sampled  phase  screen  vs.  position. 


To  employ  FFT  techniques,  we  require 

ax  -  [XZ*/N]*  ( B-3  ) 

where  N  is  the  number  of  poirts  in  the  FFT.  Then  noting  the 
periodicity  of  we  can  write  h,,,  as  a  common  sum  from  0  to 

N  -  1, 

N- 1 

hm  e-in/4  einm2/M  J  elnj2/N  ^1**  e~i2nmJ/N  (B.4) 

J»0 


which  is  in  the  :  candard  format  of  a  FFT. 
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B.3  SAMPLING  CRITERIA. 


In  order  to  accurately  calculate  the  scattered  field  by  means 
of  the  discrete  sum.  Equation  B-4,  we  have  to  adequately  sample 
the  phase  screen,  that  is,  the  phase  cannot  change  by  more  than 
n  from  one  sample  point  to  the  next, 

| <£ ( x+ax) -$( x. |  <  n,  (B-5) 


An  average  on  this  requirement  leads  to  a  bound  on  given 
by 


orMAX 


AX 


( B-6 ) 


where  the  function  F  is  defined  as 


_2  ,  l  r  dk  ♦(  1c )  ,  .  . 

F  (a>  =  J  2 n  - 2 -  ( 1-008  a  Lo*> 

a 


( B-7  ) 


For  the  power  law  spectrum,  Equation  3,  this  becomes 


F2  (a,*)  -  [l  -  [.(^|<|  (f )  K„_*(a>]/a2.  <B-8) 

The  other  necessary  condition  for  the  finite  sum.  Equation  B-4, 
to  represent  the  scattered  field  is  that  little  energy  As  scat¬ 
tered  from  the  neglected  parts  of  the  phase  screen,  that  is,  we 
can  ignore  edge  effects.  Since  the  angle  of  scatter  is 
approximately 
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(B-9) 


e 


1  dA 
k  dx' 


we  require 


1  I »(*+*.X-L._-*l*U  <1L  (B-io) 

k  ax  2  2* 

or,  in  an  averaged  sense,  a  bound  on  identical  to  Equation  B-6 . 

The  degree  to  which  our  calculations  satisfy  the  bound  Equation 
B-6  is  shown  in  Table  2  for  the  single  power  law  cases  and  in 
Table  2  for  the  two  power  law  cases.  Shown  there  is  the 
relationship  between  X^y  and  e r£,  determined  either  by  Equation 
^-9  or  by  numerically  integrating  Equation  14.  The  maximum 
Xpy  we  considered  was  =  20  with  the  corresponding  «r^  given 

in  the  Tables.  The  bound  determined  from  Equation  B-6  is  given 
in  the  Tables  as  «rMAX.  As  can  be  seen,  all  the  values  of  of 
interest  to  us  are  comfortably  small  compared  to  the  maximum 
values  so  the  conditions  on  sampling  and  edge  effects  are  well 
satisfied . 

Finally,  there  is  the  requirement  that  the  PSL  is  adequately 
sampled.  This  is  accomplished  if  the  sampled  spatial  extent, 

L,  satisfies 


L  -  5  Lq 


(B~lla) 


or 


Ak  -  1/L0. 


(B-llb) 
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TABLE  3.  SINGLE  POWER  LAW  PARAMETERS. 


s 

** 

Pi/ 

(*Ry  =  20) 

•max 

1 

11.7 

52 

264 

3/2 

30 . 8 

138 

687 

2 

81.6 

365 

2818 

5/2 

208.9 

934 

5804 

3 

501 . 8 

2244 

8208 

TABLE  4. 

TWO  POWER  LAW  PARAMETERS. 

J»F  (  in  ra) 

Ply]* 

•*  <*Ry  =  20) 

"MAX 

1386 

18 . 8 

84.3 

886 

693 

33.4 

149 . 5 

1766 

346 

72 . 6 

324 . 5 

3543 

173 

185.2 

828.4 

7150 
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This  requirement  is  met  for  all  the  single  power  law  PSD 
considerations.  For  a  few  of  the  cases  considered  for  the 
two  power  law  PSD ,  the  requirement  is  violated  by  at  most  a 
few  factors  of  two.  However,  the  very  small  k  portion  of 
the  spectrum  contributes  little  to  the  results  so  that  the 
sampling  is  adequate  even  for  these  cases.  We  did  perform  a 
few  calculations  at  larger  values  of  N  and  confirmed  that  we 
obtained  similar  results  to  those  presented. 

B.4  NUMERICAL  VALUES  OF  PARAMETERS. 

For  our  calculations,  we  assumed 

Lq  *  10  km, 

A  ^  38  10  m, 

N  =  214. 

The  probability  distributions  were  calculated  by  partitioning 
the  intensity  range  0  to  15  into  bins  of  *1  with 

Al  »  10-2. 


For  every  combination  of  PSD,  Fresnel  length  and  640 

phase  screen  realizations  were  created,  each  supplying  14384 
intensity  samples  *  Consequently,  the  total  number  of  num¬ 
bers  that  contributed  to  a  giver,  probability  density  curve 
was  about  9.2  million  and  the  minimum  probability  denrity 
measurable  is  about 


PMIN  "  10_5 * 
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Of  particular  concern  for  moderate  values  of  S4  was  the 
behavior  of  the  probability  distribution  for  small  values  of 
the  intensity.  The  behavior  here  determines  the  communica¬ 
tion  channel  behavior  for  large  values  of  SNR.  We  Investi¬ 
gated  the  small  I  behavior  in  order  to  verify  that  P(**Q) 
was  indeed  non-zero  and,  in  particular,  that  it  is  typically 
significantly  different  than  the  result  for  the  Rice  distri¬ 
bution, 

PRICE<°I  -  1_v-^-Tr'  eXP  /<  l-/T=S*l]  .  (B-12) 

1  ~S4 


Of  course,  Nakagami-m  goes  to  zero  in  this  limit.  The  pro¬ 
cedure  was  to  consider  the  interval  0  (  I  {  0.05  in  bins  of 
aI  with 


A I  *  10“4. 
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APPENDIX  C 


PROBABILITY  DISTRIBUTIONS 
FOR  SINGLE  POWER  LAW  PSD 

This  Appendix  contains  calculated  probability  distributions 
for  single  power  law  PSDs  under  conditions  where  2n  Lo  »  *F 
>>  A | .  These  distributions  are  compared  to  Nakagaml-m 
(labeled  as  N)  and  Rice  (labeled  as  R)  if  S4  is  less  than 
unity  and  to  Rayleigh  (labeled  as  SRF)  if  S4  is  greater  than 
or  equal  to  unity.  Throughout,  the  calculated  are  labeled  C. 


15 


0  i.0  2.0  3.0  H.0 


Figure  C-2„  Probability  Distribution  for 

s  “  1 .  y2  a  0.025,  S  2  =  0.095 


Figure  C-13.  Probability  Distribution  for 
s  =  1.5,  X2  =  0.01,  S  2  =  0. 


Figure  C-31.  Probability  Distribution  for 

7  7 

s  -  2,  X£  -  0.7,  *  1.12. 


[f, 


Figure  C-32o 


Probability  Distribution  for 
=  2,  X2  =  1.5,  S  2  «  1.20. 


s 


.0  i.0  2.0  3.0  4.0 

Figure  C--39,  Probability  Distribution  for 

s  =  2.5,  x2  =  0.05,  S  2  -  0.213, 
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APPENDIX  D 


PROBABILITY  DISTRIBUTIONS  FOR  SMALL 
INTENSITY  FOR  SINGLE  POWER  LAW  PSD 

This  Appendix  contains  all  the  calculated  small  intensity 
behavior  for  single  power  law  PSD.  These  distributions  are 
compared  to  Nakagaml-m  (labeled  as  N)  and  Rice  (labeled  as  R) . 
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'0.0E-0 


Figure  D- 


1.0F-02  2.0F-02  3.0E-02  4.0E-02  5.0E-02 

i 


Probability  Distributions  for  Small  I  I  >r 

2  2 

-  1,  x  -  0.1,  s/  = 


s 


0.327. 


0E-01  1.0E-02  2.0E-02  3.0E-02  4.0E-02  5.0E-02 


Figure  D- 


I 


Probability  Distribution  for  Small  I  for 
s  =  2,  X2  -  0.1,  S  2  =  0.378. 


PCI) 
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APPENDIX  E 


CALCULATED  BIT  ERROR  RATES 
FOR  SINGLE  POWER  LAW  PSD 


This  Appendix  contains  all  the  calculated  bit  error  rates  for 
single  power  law  PSD.  The  six  modems  considered  are:  CPSK, 
APSK,  DBPSK,  BFSK ,  QFSK  and  8-ARYFSK .  The  curves  are  labeled 
as  follows: 


AWGN 

a 

b 

c 

d 

e 

f 

SRF 


Additive  white  Gaussian  noise 


0.01 


=  0.025 


=  0.05 


=  0.1 

=  0.25 


=  0.4 


Slow  Rayleigh  Fading 
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Figure  E-5.  BER  for  QFSK 


Figure  E-26.  BER  for  APSK  s 
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APPENDIX  F 


CALCULATED  BIT  ERROR  RATES  FOR  DBPSK 
COMPARED  TO  NAKAGAMI-M  AND  RICE 
FOR  SINGLE  POWER  LAW  PSD 


This  Appendix  contains  the  calculated  bit  error  rate  for  DBPSK 
compared  to  Nakagami-m  (labeled  as  N)  and  Rice  (labeled  as  R) 
for  single  power  law  PSD.  Throughout,  the  calculated  results  are 
labeled  c  . 
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APPENDIX  G 


PROBABILITY  DISTRIBUTIONS 
FOR  TWO  POWER  LAW  PSD 

This  Appendix  contains  all  the  calculated  probability  distri¬ 
butions  for  two  power  law  PSD.  These  distributions  are  com¬ 
pared  to  Nakagami-m  (labeled  as  N)  and  Rice  (labeled  as  R)  if 
S4  is  less  than  unity  and  to  Rayleigh  (labeled  as  SRF)  if  is 
greater  than  or  equal  to  unity.  Throughout,  the  calculated 
results  are  labeled  C. 
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APPENDIX  H 


PROBABILITY  DISTRIBUTIONS  FOR  SMALL 
INTENSITY  FOR  TWO  POWER  LAW  PSD 


This  Appendix  contains  all  the  calculated  small  intensity 
behavior  for  two  power  law  PSD.  These  distributions  are 
compared  to  Nakagaml-m  (labeled  as  N)  and  Rice  (labeled  as  R) . 
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Probability 

a  «  2,  X2  - 


Distribution  fo. 

0.4,  S42  *  0.85 
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APPENDIX  J 


CALCULATED  BIT  ERROR  RATES  FOR  DBPSK 
COMPARED  TO  NAXAGAMI-M  AND  RICE 
FOR  TWO  POWER  LAW  PSD 

This  Appendix  contains  the  calculated  bit  error  rate  for  DBPSK 
compared  to  Nakagami--m  (labeled  as  N)  and  Rice  (labeled  as  R) 
for  two  power  law  PSD.  Throughout,  the  calculated  data  is 
labeled  C . 
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SPACE  &  NAVAL  WARFARE  SYSTEMS  CMD 
ATTN:  CODE  3101  T  HUGHES 
ATTN:  PD50TD 
ATTN:  PD50TD1  G  BRUNHART 
ATTN:  PME  106-4  S  KEARNEY 
ATTN:  PME-106  FWDIEDERICH 

THEATER  NUCLEAR  WARFARE  PROGRAM  OFC 
ATTN:  PMS-42331F(D  SMITH) 


SANDIA  NATIONAL  LABORATORIES 
ATTN.  D  HARTLEY 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  ADTHORNBROUGH 
ATTN-  CODE  9014  RBACKSTROM 
ATTN.  D  DAHLGREN 
ATTN:  ORG  1231  TP  WRIGHT 
ATTN:  ORG  9114  WD  BROWN 
ATTN:  SPACE  PROJECT  D!V 
ATTN:  TECH  LIB  3141 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  OSWR/NED 
ATTN:  OSWR/SSD  FOR  L  BERG 


DEPARTMENT  OF  THE  AIR  FORCE 

AIR  FORCE  CTR  FOR  STUDIES  &  ANALYSIS 
ATTN:  AFCSA/SAMIR  GRIFFIN 
ATTN:  AFCSA/SASC 

AIR  FORCE  ELECTRONIC  WARFARE  CENTER 
ATTN.  LT  M  MCNLELY 

AIR  FORCE  SPACE  SYSTEMS  DIVISION 
ATTN:  YA 
2 CYS  ATTN:  YN 

AIR  FORCE  TECHNICAL  APPLICATIONS  CTR 
ATTN:  TN 

AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL-LSE 

HQ  AWS,  DET  3  (CSTC/WE) 

ATTN:  WE 

SECRETARY  Or  AF/AQQS 
ATTN:  AF/RDQI 

STRATEGIC  AIR  COMMAND/XRFS 
ATTN.  XRF$ 

WEAPONS  LABORATORY 
ATTN:  NTCA 


DEPARTMENT  OF  COMMERCE 
ATTN:  E  MORRISON 
ATTN:  J  HOFFMEYER 
ATTN:  WUTLAUT 

U  S  DEPARTMENT  OF  STATE 

ATTN:  PM/TMP  RM  7428 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN:  A  EIGHTY 
ATTN:  BPPURCFLl 
ATTN:  C  CREWS 
ATTN:  C  RICE 
ATTN:  G  LIGHT 
ATTN:  M  ROLENZ 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
ATTN  SECURITY 

ATLANTIC  RESEARCH  SERVICES  CORP 
ATTN:  R  MCMILLAN 

ATMOSPHERIC  AND  ENVIRONMENTAL  RESEARCH  INC 
ATTN:  MKO 

AUSTIN  RESEARCH  ASSOCIATES 
ATTN:  J  THOMPSON 
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UNA-Ttm-ISQ  (DL  CONTINUED) 


AUTOMETRIC  INCORPORATED 
ATTN  C LUCAS 

BDM  INTERNATIONAL  INC 
ATTN:  L JACOBS 

BERKELEY  RSCH  ASSOCIATES,  INC 
ATTN:  C  PRETTIE 
ATTN  J  WORKMAN 
ATTN:  S  BRECHT 

BOEING  CO 

ATTN:  GHALL 

CALIFORNIA  RESEARCH  &  TECHNOLOGY,  INC 
ATTN:  M  ROSENBLATT 

CHARLES  STARK  DRAPER  LAB.  INC 
ATTN:  ATETEWSK! 

COMMUNICATIONS  SATELLITE  CORP 
ATTN:  G  HYDE 

CORNELL  UNIVERSITY 

ATTN.  D  FARLEY  JR 
ATTN  M  KELLY 

ELF.CTROSPACE  SYSTEMS,  INC 
ATTN:  P  PHILLIPS 

EOS  TECHNOLOGIES,  INC 
ATTN:  B GABBARD 
ATTN:  W  LELEVIER 

GENERAL  RESEARCH  CORP  INC 
ATTN:  J  EOLL 

GEO  CENTERS,  INC 

ATTN:  E  MARRAM 

GRUMMAN  AEROSPACE  CORP 
ATTN:  J  DIGLlO 

GTE  GOVERNMENT  SYSTEMS  CORPORATION 
ATTN:  W  I  THOMPSON,  III 

HARRIS  CORPORATION 
ATTN:  E  KNICK 

HSS,  INC 

ATTN:  D HANSEN 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  E BAUER 
ATTN:  H  WOLFHARD 

JS  LEE  ASSOCIATES  INC 
AT  iN  DR  J  LEE 

JAYCOR 

ATTN:  J  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
ATTN.  CMENG 
ATTN:  JD  PHILLIPS 
ATTN:  R STOKES 
ATTN:  T  EVANS 


KAMAN  SCIENCES  CORP 
ATTN:  E  CONRAD 
ATTN;  GDiTTBERNER 

KAMAN  SCIENCES  CORPORATION 
ATTN:  BGAMBILL 
ATTN:  DASIAC 
ATTN:  R RUTHERFORD 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  J  HENL  EY 
ATTN:  J  KUMER 
ATTN:  R SEARS 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  DKREJCI,  YOOl 

LTV  AEROSPACE  &  DEFENSE  COMPANY 
2  CYS  ATTN:  LIBRARY 

M  I  T  LINCOLN  LAB 

ATTN:  D  TOWLE 
ATTN:  I  KUPIEC 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN.  H  VON  STRUVE  III 

MAXIM  TECHNOLOGIES,  INC 
ATTN:  BRIDGEWAY 
ATTN:  JSO 

MCDONNELL  DOUGLAS  CORP 
ATTN:  TCRANOR 

MCDONNELL  DOUGLAS  CORPORATION 
ATTN:  J GROSSMAN 
ATTN:  RHALPRIN 

METATECH  CORPORATION 
ATTN:  R SCHAEFER 
ATTN:  WRADASKY 

METEOR  COMMUNICATIONS  CORP 
ATTN:  R LEADER 

MISSION  RESEARCH  CORP 
ATTN:  IKENNEALY 
ATTN:  R  ARMSTRONG 
ATTN:  R  LARKIN 
ATTN:  S  BELANGER 
ATTN:  W  WHITE 

MISSION  RESEARCH  CORP 
ATTN:  BR  MILNER 
ATTN:  D  ARCHER 
ATTN:  DKNEPP 
ATTN:  D  LANDMAN 
ATTN:  FFAJEN 
ATTN:  F  GUIGLIANO 
ATTN:  GMCCARTOR 
ATTN:  KCOSNER 
ATTN:  M  FIRESTONE 
ATTN:  RBIGONI 
ATTN:  R  BOGUSCH 
ATTN:  R  DANA 
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ONA-TR-89  150  (DL  CONTINUES)) 


ATTN:  R  HENDRICK 
ATTN:  RKILB 
ATTN:  SGUTSCHE 
ATTN:  TECH  INFO  CENTER 
ATTN:  TECH  LIBRARY 

MITRE  CORPORATION 

ATTN:  M  R  DRESP 

MITRE  CORPORATION 

ATTN:  J  WHEELER 
ATTN:  MHORROCKS 
ATTN:  RCPESCI 
ATTN:  W  FOSTER 

NORTHWEST  RESEARCH  ASSOC,  INC 
ATTN:  EFREMOUW 


RJO  ENTERPRISES/POET  FAC 
ATTN:  A  ALEXANDER 
ATTN:  W BURNS 

S-CUBED 

ATTN:  C NEEDHAM 
ATTN:  T  CARNEY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  C  SMITH 
ATTN:  D  HAMLIN 
ATTN.  D SACHS 
ATTN:  ESTRAKER 
ATTN:  LLINSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J COCKAYNE 


PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  E  FIELD  JR 
ATTN:  F  THOMAS 
ATTN:  H  BRODE 

PHOTOMETRICS,  INC 

ATTN:  I  L  KOFSKY 

PHYSICAL  RESEARCH  INC 
ATTN:  W.  SHIH 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DTELAGE 
ATTN:  M  CROSS 

SRI  INTERNATIONAL 

ATTN:  WCHESNUT 
ATTN:  WJAYE 

STEWART  RADIANCE  LABORATORY 
ATTN:  RHUPPI 


PHYSICAL  RESEARCH  INC 
ATTN:  HFITZ 
ATTN:  PLUNN 

PHYSICAL  RESEARCH,  INC 
ATTN:  RDELI3ERIS 
ATTN:  T STEPHENS 


TELECOMMUNICATION  SCIENCE  ASSOCIATES 
ATTN:  R BUCKNER 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  J  WOLFSBERGER,  JR 

TOY0N  RESEARCH  CORP 
ATTN:  J  ISE 


PHYSICAL  PESEARCH,  INC 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  WSCHLUETER 


TRW  INC 

ATTN:  DRDGRYBOS 

ATTN:  R  PLEBUCH/HARDNESS  &  SURV  LAB 
ATTN:  H  CULVER 


R  &  D  ASSOCIATES 


ATTN: 

C  GRtlFINGER 

ATTN. 

F  GILMORE 

ATTN: 

G  HOYT 

2  CYS 

ATTN: 

L  L DERAAD 

ATTN 

M GANTSWEG 

2  CYS 

ATTN: 

M  K  GROVER 

R  &  0  ASSOCIATES 

ATTN:  J  WALTON 


TRW  SPACE  &  DEFENSE  SYSTEMS 
ATTN:  D  V  LAYTON 

USER  SYSTEMS,  INC 

ATTN:  S  W  MCCANDLESS  'R 


UTAH  STATE  UNIVERSITY 

ATTN:  K  BAKER,  DIR  ATMOS  &  SPACE  SCI 
ATTN:  L  JENSEN,  ELEC  ENG  DEPT 


RAND  CORP 

ATTN:  C  CRAIN 
ATTN:  E  BEDROZIAN 


VISIDYNE,  INC 

ATTN:  J CARPENTER 

FOREIGN 


RAND  CORP 

ATTN:  B  BENNETT 


FOA  2 

ATTN.  BSJOHOLM 


FOA  3 

ATTN:  T  KARLSSCN 
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